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automatic  control  systems  for  the  vertical  motion  of  a 
fictitious  submarine.  A  linearized  mathematical  model 
is  derived  from  a  non-linear  model  in  six  degrees  of 
freedom.  Both  designs  are  simulated  and  evaluated  with 
respect  to  performance,  simplicity  of  design 
procedure,  and  grade  of  complexity.  The  controller, 
designed  via  the  steady-state  decoupling  method,  is 
implemented  in  the  non-linear  model  and  tested  under 
various  operating  conditions. 
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Characteristic  for  a  submarine  is  translation  in 
dimensions,  where  depth  control  is  of  essential  impor 
The  automatic  control  of  motion  in  the  vertical  plane 
depth  and  pitch  control,  is  the  subject  of  this  thesis 
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This  thesis  investigates  two  possible  design  methods, 
steady-state  and  total  decoupling.  As  a  basis  for  both 
designs  a  linear  model  is  developed  and  validated. 

Steady-state  decoupling  is  achieved  by  designing  a 
cascaded  diagonal  compensator  matrix,  using  classical 
sfngle-loop  techniques.  Based  on  the  desired  response 
characteristics  the  total  decoupling  scheme  leads  to  a 
cascaded  compensator  matrix  with  off-diagcnal  terms.  Both 
designs  are  simulated  and  evaluated  with  respect  to 
performance,  simplicity  of  design  procedure,  and  grade  of 
complexity. 
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II.  SIMOLATION  H0D2LS 


A.  STANEARC  MODEL 

A  set  of  "  Standard  Equations  of  Motion  for  Submarine 
Simulation  "  was  developed  by  NSRDC  (Ref.  1).  Equations  and 
additional  auxiliary  equations  are  repeated  in  Appendix  A. 

These  equations  are  referred  to  a  right-hand  orthogonal 
system  of  moving  axes,  fixed  in  the  body  with  its  origin 
looated  at  the  oenter  of  mass  (CG)  of  the  body.  The 
x-z-plane  is  the  principle  plane  of  symmetry  (  vertical 
center  plane  for  submarines  ) ;  the  x-axis  is  parallel  to  the 
baseline  of  the  body.  Angular  velocity  components,  forces, 
and  moments  are  shown  in  Fig.  1. 

The  equations  are  written  in  a  form  utilizing 
nondimensionalized  hydrodynamio  coefficients  and  are 
applicable  to  the  rigid  body  motions  of  submarines  and  other 
submerged  vehicles.  Complete  sets  of  hydrodynamic 
coefficients  have  been  determined  by  NSRDC.  The 
coefficients  for  a  fictitious  submarine  used  in  this  thesis 
stem  from  a  submarine  simulation  program  developed  in 
Ref.  2.  A  notation  for  the  various  terms  in  the  equations 
including  the  hydrodynamic  coefficients  is  given  in  Appendix 

B.  Numerical  values  of  coefficients  and  parameters  are  given 
in  Appendix  C.  All  coefficients  of  the  terms  containing 
(  -1)  and  OC  are  set  to  zero.  These  terms  reflect  the 
incremental  changes  in  forces  and  moments  -  generated  by 
propeller  rpo  -  due  to  either  over  or  underpropulsion .  For 
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Figure  1  -  SOBMARINE  AXES 
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the  moderate  changes  in  ahead  speed  involved  in  most  normal 
maneuvers  all  terms  can  be  neglected.  The  coefficients  used 
in  this  thesis  apply  only  to  the  deeply  submerged  case,  free 
of  near- surf ace ,  bottom,  and  wall  effects. 

As  Bef.  1  points  out,  correlation  studies  and 
preliminary  comparisons  have  shown  that  the  "Standard 
Equations"  together  with  a  given  set  of  coefficients  yield 
accurate  predictions  of  normal  maneuvers  in  submerged  ahead 
motion.  The  "Standard  Equations  of  Motion"  will  be  referred 
to  as  "Noc-iinear  Model". 


In  Bef.  2  H 
Digital  Simulatio 
equations  was  so 
Cramer's  rule  dev 
cofactors  fcr  th 
this  thesis  here 
parameters.  The  v 
shewn  by  Bef.  2. 
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Ived  simultaneously  by  an  approach  using 
eloped  and  programmed  by  fief.  3.  The 
e  set  of  hydrodynamic  coefficients  used  in 
calculated  by  Drurey  and  are  used  as 
alidity  of  the  DSL  implementation  has  been 
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B.  DERIVATION  OF  THE  LINEARIZED  MODEL 


This  section  describes  the  assumptions/  procedures,  and 
tests  to  validate  the  linear  model. 


1  •  Assunrtions 


♦  Forward  speed  is  constant 


An  attempt  to  linearize  about  the  axial  speed,  u, 
which  affects  nearly  every  term  in  the  standard  equations, 
would  imply  a  great  deal  of  complexity.  Thus  the  forward 
speed  was  assumed  to  be  constant.  This  means  that  either  the 
forward  speed  variations  are  small  and  can  be  neglected  or 
if  variations  are  large  they  impose  a  control  problem  and 
will  be  solved  accordingly. 

*  Boll  angle  is  small 


In  submarine  maneuvering  large  roll  angles  occur 
only  in  the  fcrm  of  "  Snap  Roll  ",  i.e.  high  speed  plus  hard 
over  rudder.  Under  normal  conditions  depth-changing  and 
-keeping  maneuvers  produce  very  small  roll  angles.  Therefore 
the  roll  angle  can  be  neglected. 

*  Cross  products  of  inertia  may  be  neglected 


I  his 

simulations 
submarines  is 


assumption  is  common 
because  the  hull  and 
approximately  symmetric. 


to  all  submarine 
interior  layout  of 
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*  All  terms  involving  W  are  neglected 

i 

This  is  possible  because  the  submarine  was  assumed 
to  be  in  triii. 

♦  The  submarine  can  be  approximated  by  a  linear  model 

The  linear  equations  are  assumed  to  be  a  valid 
representation  in  the  depth- keep ing  and  depth-changing  mode 
of  operaticn  respectively.  This  assumption  is  not  obvious 
but  will  be  verified  via  simulation. 

These  assumptions  lead  to  the  decoupling  of  the  vertical 
plane  from  the  remaining  equations. 

2*  estimation  of  state  equations 

As  this  thesis  is  only  concerned  with  the  vertical 
plane,  the  linearization  of  the  horizontal  plane  is  omitted. 

Thus  the  linear  equations  are  (FORTRAN-Notation,see 
Appendix  C)  : 

4  3 


^♦WDCT  -  DC*il*Q 

=  L  ♦ZCDOT’^QDOT  + 

L  ♦ZWDOT+WDCT  + 

3 

L  *ZQDOT*UC*Q  + 

2 

L  ♦Zy*UC*W 

2  2 

L  *00  ♦  (ZDS*DS 

f  ZDB*DB) 

(1) 

IY*CDOT 

5 

=  L  *aQD0T*QD0T  + 

4 

L  ♦.aQ*UC*Q 

4 

L  *MWDOT*MDOT  + 

3 

L  ♦MW*UC*W 

+ 

3  2 

L  *UC  (MDS*DS  + 

MDB*DB)  + 

E*ZE*PITCH  (2) 


13 


m 

r 

t  ^  1 


I 


\ 


The  linearized  auxiliary  equations  are  : 


IPIDCT 


LZEOT 


=  LQ 

=  W  -  UC=*‘LPITCH 


(3) 

(4) 


Solving  equation  (1)  and  (2)  for  WDOT  and  LQDOT 
respectively  and  substituting  into  each  other  yields  the 
following  relations  ; 


'hdot  “ 

a 

11 

a  a 

12  13 

W 

~b 

1  1 

b 

12 

ICDOT 

3 

a 

21 

a  a 

22  23 

LQ 

b 

21 

b 

22 

1C 

a  ^ 

<4 

a  a 

1  3  2  33 

LPITCH 

b 

31 

b 

32 

DS 

DB 


(5) 


which  can  he  written  as 


where 


1 1 


12 


13 


21 


X  =  Ax  +  BR 

-irf  ^  *\j 


2  1  5 
[L  *ZW*UC  +  L  *ZQDOT*aW*UC/ (lY-L  *flQDOT)]/^ 

3  3  5 

[L  *ZC*UC  +  K*UC  +  L  ♦ZQDOT=^MQ*UC/ (lY-L  *MQDOT)  ]/^ 

4  5 

■[L  =*ZCEOI*ZB*B/ (lY-L  *MQDOT)]/j^ 

3  6  3 

[L  *tlS’*UC  +  L  ♦MWDOT*ZW*UC/ (M-L  *ZWDOT)  ]/£ 


(6) 


4  7  4 

a  _  =  [L  ^I1Q=*UC  +  L  *MWDOT*ZQ*OC  +  L  *MHDOT»M*UC/ 
2  2 


(M-L  +2KD0T)  ]/£ 

I  =-B*ZbA 
23  ^ 


14 


• 


.1 


t 


i 

f 


I 


ft 


fit 


# 


n 

1 
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a  =0.0 

31 

a  =  1.0 

32 

a  =0.0 

33 


2  2  7  2  5 

=  [L  *ZDS*UC  +  L  ♦ZQDOT*MDS*UC  /(lY-L  *flCDOT)  ]/ 

2  2  7  2  5 

=  [L  ♦ZEE*UC  +  L  *ZQDOT*MDB*UC  /(lY-L  ♦MCDOT)  ]/ 

3  2  6  2  3 

b^^  =  [L  ♦MDS^UC  +  L  *MWDOT*ZDS*UC  /(M-L  *ZWDOT)  ]/ 

3  2  6  2  3 

b  =  [L  =*MDE*OC  +  L  *MWDOT*2DB#UC  /  (M-L  ♦ZWDOT)  ]/ 
22 

b  =0.0 

31 

b  =0.0 

32 

with 


3  8  5 

i^=  a  -  L  *ZHDOT  -  L  *ZQDOT*MWDOT/ (I  Y-L  ^iJlCDOT) 

5  8  3 

£=  lY  -  L  ♦MQDOT  -  L  ♦MWDOT*ZCDOT/ (H  -  L  *ZWDOT) 


Vector  equation  (5)  describes  the  state  variable 
representaticn  of  the  linearized,  vertical  plane  equations 
of  motion.  However  the  state  vector  does  not  contain  ZDOT 
but  W.  W  represents  the  component  of  u  in  the  z-direction, 
in  order  to  make  the  depth  a  state  variable  one  wants  to  use 
LZDCT  which  represents  the  depth  chanqe  in  the  z-direction. 
Thus  the  linearized  auxiliary  equation  (4)  is  used  for  the 
transformation 


LZDCT  =  tf  -  UC»LPITCH 
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By  definition 


LZDCT 

w 

= 

LQ 

IX 

II 

LQ 

LPIICH 

LPITCH 

Using  the  linear 
variables  is  formed 


operator  M, 


the 


new  set  of 


state 


/N 


X 


HX 


where 


a 


1  0  -uc 

0  1  0 
0  0  1 


equations 


(6)  become 


X  +  BR 

then 


-1 

M  X  =  AM 


X 

=  MAM  X  +  MBR 

but 

/V 

^  A/ 

✓N 

MB 

=  B  =  B 

thus 

-j 

A/ 

X  =  AkH  X  +  3R 


which  can  he  rewritten  as 


X  =  Ax  +  BR 


(7) 
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^  A  i 


’  •  a 


1  ^ 


i 


where 


A 

A  = 


1 1 

12 

13 

a 

a 

[a 

2  1 

22 

21 

_0. 0 

1 . 0 

te 

written 

as 

23 


0.0 


A 

A  = 


2  n 

A11+0C  A12+UC  [A13+A1U*UC  ] 

2 

A21+CC  A22*UC  [ A23*UC  +A24  ] 


and 


B 

rxf 


=  QC 


.0 

1 . 

0 

o 

• 

o 

Ell 

B12 

E21 

E22 

C.C 

o.<^ 

Using 

the  numerical  values 

for  the 

hydrodynamic 

coefficients 

given  in  Appendix  C 

and  letting 

UC  still  be 

variable,  the 

constants  of  the  A  and 

B  matrices 

ar  e 

A1  1 

=  -1.72£ 

-3 

10 

A12 

=  -0.706 

2 

A13 

=  C.0128S 

A14 

=  -  1.728 

-3 

10 

A2  1 

=  1.684 

-5 

10 

A22 

=  -6.365 

-3 

10 

A23 

=  1.884 

-5 

10 

A24 

=  -2.522 

-3 

10 
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BIT 

=  -6.666 

-4 

10 

B12 

=  -3.871 

-4 

10 

B2  1 

=  -  1.465 

-5 

10 

E22 

=  3.19C 

-6 

10 

C.  VALICATICN  OF  LINEAR  MODEL 

The  ctjective  of  this  section  is  to  compare  the  dynamics 
of  the  standard  model  with  the  developed  linear  model. 

In  order  to  compare  both  models,  they  should  be  in  the 
same  initial  state,  i.e.  besides  the  identical  initial 
corditions  tcth  models  have  to  be  ”  in  trim  One 
linearizing  assumption  was  that  the  model  will  be  in  trim  at 
all  times.  Therefore  the  standard  model  was  brought  into 
"neutral  conditions"  for  the  given  test  speeds,  in  order  to 
avoid  additional  forcing  terms.  This  procedure  will  be 
described  explicitly  in  a  later  chapter. 

Any  total  system  response  can  be  viewed  as  having  two 
components:  forced  and  free.  The  forced  component  includes 
the  complex  frequencies  of  the  forcing  functions,  the  free 
component  the  complex  or  natural  frequencies  of  the  system. 
Both  response  components  were  checked  via  simulation. 


18 


J 


1 


Initial  condition  response 

It  has  expected  that  for  small  perturbations  the 
deviations  between  both  models  should  be  small.  Therefore 
initial  conditions  of  2.5°  in  pitch  where  tested  first. 
The  initial  value  of  the  second  variable  LZDOT  was  defined 
as  -UC*sin  (initia 1  pitch),  based  on  the  auxiliary  equation 
for  depth  change. 

lest  runs  over  400  sec  in  the  speed  range  3-15 
knots  were  performed  simultaneously  for  both  models. 
Table  01  shows  the  max.  differences  between  the  model 
responses.  Fig.  2-11  are  a  representative  selection  of  the 
depth  change  and  the  pitch  behavior.  In  each  figure  curve  1 
represents  the  non-linear,  curve  2  the  linear  model. 
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Tdble  01  -  Initial  condition  response  to  2.5  °  pitch 
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As  one  can  readily  observe,  all  deviations  are  small 
for  this  initial  condition,  which  indicates  that  the 
dynamics  are  nearly  identical  for  small  perturbations. 

Ihe  maximum  pitch  angle  expected  in  normal 
operations  is  limited  to  about  45  Therefore  this  initial 
condition  was  simulated  as  well.  As  for  the  initial 
condition  of  2.5°  the  max.  deviations  are  summarized  in 
Table  02.  In  Fig.  12-23  curve  1  represents  the  non-linear, 
curve  2  the  linear  model. 


Ihe  second  set  of  initial  conditions,  which 
represents  rather  large  perturbations,  still  leads  to  very 
similar  dynamic  behavior,  but  the  deviations  increased 
substantially.  This  was  expected  as  the  angle  approximation 

sin  0=0 


is  rather  crude  for  angles  of  this  magnitude.  In  addition 
to  that  the  constant,  speed  assumption  is  not  as  valid  as  for 


IC  =  2.5°. 
observed 


For  both  sets  of  initial 
that  increasing  speed  tends 


condition  it  is 
to  increase  the 


deviations  between  both  trajectories. 
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Table  02  -  Initial  condition  response  to  45.  °  pitch 
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2*  response 


As  the  models  considered  are  deeply  submerged  and  in 
trim,  the  planes  exert  the  only  relevant  forces.  It  was 
considered  best  to  test  the  response  of  both  models 
(excluding  actuators)  to  sinusoidal  forcing  functions, 
because  in  the  controlled  system  the  rudder  excursions  due 
to  depth  and  pitch  changes  are  expected  to  be  sinusoidal. 

The  following  forcing  functions  were  chosen  for  both 
planes  separately  and  three  different  speeds  (3,  9,  15 
kncts)  ; 
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Table  03  -  Speed 
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Table  04  -  Speed 
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*  3  knots 

The  response  to  the  sternplane  showed  an  incr 
deviation  in  pitch  and  depth  change  with  incr 
amflitude.  With  decreasing  frequency  the  gai 
increased.  The  percent  deviation  stayed  about  the  sane 

The  response  to  the  fairwater  plane  shows  th 
behavior  except  that  the  gain  ratio  DS/DB  =  3.5  for 
and  1.6  for  depth  change.  The  absolute  oagnitude 
deviations  is  considered  negligibly  small. 

This  can  te  expected,  because  the  pitch  involv 
relatively  sttall  and  the  speed  is  about  constant 
decrease  in  speed  is  .38  kn  for  35®  plane  deflection) . 

*  S  knots 

The  speed  decreases  substantially  (up  to  2.2  k 
large  place  angles,  which  causes  bigger  deviations  in 
and  depth  change  (up  to  12.7°). 

But  it  cculd  be  observed  that  the  linearized  model  was 
a  good  approximation  for  small  perturbations. 

*  15  knots 

At  this  speed  the  decrease  of  a  was  even  bigg 
to  32%  at  35®  plane  deflection) ,  whicn  resulted  in 
deviations  in  both  variables.  For  small  perrurbatio 
approximation  is  still  very  good. 

Thus  it  has  been  shown  that: 

-  Approximation  by  linear  model  is  valid  for 
perturbations  at  all  speeds. 

-  The  dynamics  of  both  models  are  compatible 

In  addition  it  has  been  observed,  tha 
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linearized  acdel  is  still  valid  for  large  perturbations 
applied  ever  a  short  period  of  time  (initial  condition 
response) .  When  large  plane  angles  were  applied  over  a  long 
period,  the  speed  of  the  standard  model  decreased 
substantially,  which  made  the  constant  speed  assumption 
invalid . 


3.  Actuator  model 
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which  has  the  Laplace  transform  of 

1  1 
s  (s  +  oL) 

Thus  the  actuator  dynamics  will  be  represented  by: 
1 

(s+d) 
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ot,  represents  the  inverse  of  the  time  constant, 
which  has  to  be  determined.  As  the  rate  of  the  rudder 
deflection  depends  on  the  magnitude  of  the  desired  angle  the 
maximum  rate  limit  is  determined  by  the  mechanical  limits  of 
the  rudder,  which  are  assumed  to  be  ±  35  °  for  this  model. 
It  was  assumed  that  the  maximum  deflection  is  reached  in  5 
time  constants.  As  the  needed  time  is  about  3  to  3.5  sec, 
the  time  constant  was  determined  to  be  2/3  sec. 


The  actuator  time  constant  is  fast  compared  to  those 
of  the  system.  Therefore  the  actuators  will  be  neglected  in 
the  design  of  the  compensators,  but  will  be  included  in  the 
simulation  of  the  compensated  system. 
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III.  DESIGN  SPECIFICATIONS 


The  objective  is  to  design  an  automatic  controller  for 
maneuvering  in  the  vertical  plane,  i.e.  depth  and  pitch 
control.  The  submarine  is  considered  deeply  submerged  and 
therefore  free  from  all  external  disturbances  like  surface-, 
bottom-,  and  wall-effects.  The  dynamics  for  a  depth  change 
include  the  depth  and  pitch  dynamics.  The  deptn  response  to 
an  ordered  depth  change  should  basicly  look  like  a  second 
order  overdamped  or  highly  damped  system  response.  The  pitch 
shall  be  kept  as  close  as  possible  to  the  ordered  pitch, 
which  under  regular  conditions  will  be  zero.  This  submarine 
has  no  inertial  guidance  system  (as  for  example  the  German 
coastal  submarines) .  Sate  information  is  nor  available  and 
the  only  states  to  be  used  as  feedback  are  depth,  pitch,  and 
speed.  Cue  to  this  limited  instrumentation  the  controller 
will  have  to  use  cascaded  filters. 


*  Time  requirements. 

Depth  changes  at  even  keel  for  speeds 


-  <  9  kn 

10 

f t  in  < 

120  sec 

100 

ft  in  < 

240  sec 

-  >  9  kn 

10 

ft 

in  < 

60  sec 

100 

ft 

in  < 

100  sec 
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Overshoot  requirements 


For  depth  changes 

<  1Q0  ft  <  5  % 

>  100  ft  <  5  ft 

♦  Pitch  deviations 

Always  less  than  ±  2  °  from  ordered  pitch. 

*  Mechanical  constraints 

-  plane  angle  limits  according  to  diagram 


Speed  in  kn 


maximum  possible  pitch  <  45  ° 
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IV. 


DESIGN 


A.  TRANSPER  FUNCTION 


The  design  approach  chosen  for  this  thesis  requires  the 
system  description  by  the  transfer  function  matrix.  Matrix 
equation  (7)  is  the  state  variable  representation  cf  the 
system.  The  signal  flowgraph  equivalent  for  this  system  is 
given  below.  An  additional  integration  is  included,  because 
the  depth  is  a  desired  output. 
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The  input-output  relations  of  interest  are  DEPTH/DS, 
DEETH/DB,  PITCH/DS ,  and  PITCH/DB,  because  depth  and  pitch 
will  be  the  measurable  quantities  and  will  be  used  for  the 
control . 

Applying  Mason's  gain  rule  the  input-curput  relations 
have  the  form 

2 

Y(s)  N2*s  +  N1*s  +  NO 

=  0  3  2 

X  (s)  EU  +  s  +  D3*s  +  D2*s  +  D1*s  *  DO  (8) 


The  coefficients  for  equation  (8)  are; 


ZECT/DS 

• 

N2 

= 

2 

B  1 1*UC 

N1 

= 

3 

(  A22*B  1  1  -  A  12»B2  1)  *00 

NO 

= 

2  2  2 
(E21*(A13  +  A14  +  UC  )  -  B11=*(A23*0C  +  A24))*UC 

D4 

= 

1 

D3 

= 

-(All  +  A22) *UC 

D2 

= 

2  2  2 
A11^A22*UC  -  (A23*UC  +  A24)  -  A12»A21*DC 

D1 

= 

2  2 
A11*(A23<‘0C  +  A24)»UC  -  A21»UC*(A13  +  A14*UC  ) 

DO 

= 

o 

• 

o 

ZDOT/DB 

: 

N2 

2 

B  12*0C 

N1 

= 

3 

-(A22*B12  -  A12*322)*UC 

NO 

= 

2  2  2 
(E22*(A13  +  A14+0C  )  -  B12*(A23*UC  +  A24))*UC 

D4 

= 

1 

D3 

= 

-(All  +  A22)*UC 
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D2  = 

2  2  2 
&11*A22*UC  -  (A23»UC  +  A24)  -  A12*A21»UC 

D1  = 

2  2 
A11*(A23*UC  +  A24) *UC  -  A21*UC*(A13  +  A14*aC  ) 

DO  = 

o 

• 

o 

PIICH/DS 

• 

• 

N1  = 

2 

E21+UC 

HO  = 

3 

(A2  1*B1  1  -  A1  1»32  1)  ^UC 

D3  = 

1 

D2  = 

-(All  +  A22)  *UC 

D1  = 

2  2  2 
A11*A22*UC  -  (A23*UC  +  A24)  -  A12*A21*UC 

DO  = 

2  2 
A11*(A23*UC  +  A24) *00  -  A21^UC*(A13  ♦  A14*UC  ) 

PITCH/DB 

• 

N1  = 

2 

E22*UC 

NO  = 

3 

(a21*B12  -A1  1*B22) *00 

D3  = 

1 

D2  = 

-(All  +  A22)*UC 

D1  = 

2  2  2 
A11*A22*UC  -  (A23+UC  +  A24)  -  A12*A21*0C 

DO  = 

2  2 
A11»(A23*0C  +  A24)  *UC  -  A21*0C*(A13  +  A14^(3C  ) 

Table 

06  through  10  show  the  root-locations  for  the 

speed  range  1  kn  -  15  kn.  The  characteristic  roots  move 
closer  to  the  origin  and  the  damping  reduces  with  decreasing 
speed.  This  indicates  that  the  overall  response  of  the 
uncompensated  submarine  will  be  much  slower  for  3  kn  than 
for  9  kn  cr  15  kn. 
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Table  06  -  Characteristic  roots 


Speed 

Root 

Complex 

roots 

1 

-0.2774 

-2 

10 

-0.5452 

10  ^ 

) 

± 

j 

0.4978 

-1 

10 

2 

-0.5650 

-2 

10 

-0. 1085 

-1 

10 

1 

± 

j 

0.484  1 

-1 

10 

3 

-0.876  e 

-2 

10 

-0.1613 

-1 

10 

1 

± 

j 

0.  4604 

-1 

10 

4 

-0.1235 

-  1 

10 

-0.2114 

- 1 

10 

1 

± 

j 

094245 

-1 

10 

5 

-0. 1705 

-  1 

10 

-0.2563 

-1 

10 

1 

± 

j 

0.3716 

-1 

10 

6 

-0.2527 

-  1 

10 

-0.2835 

-1 

10 

1 

± 

j 

0.2908 

-1 

10 

7 

-0.4872 

-1 

10 

-0.2347 

-1 

10 

1 

± 

j 

0.2116 

-1 

10 

8 

-0.7068 

-  1 

10 

-0. 1933 

-1 

10 

1 

± 

j 

0. 2032 

-1 

10 

9 

-0.88  1  2 

- 1 

10 

-0. 1744 

-1 

10  ± 

j 

0.2014 

-1 

10 

10 

-0. 1038 

-0.1643 

-1 

10  1 

j 

0.  1999 

-1 

10 

1  1 

-0.  1  186 

-0.1 586 

1  0 

I 

± 

j 

0.1982 

-1 

10 

12 

-0. 1 328 

-0. 1557 

-1 

10  ± 

n 

J 

0.1963 

-1 

10 

13 

-0  .  1467 

-0. 1 546 

10 

1 

± 

3 

0.1940 

-1 

10 

14 

-0.1603 

-0.  1549 

10 

1 

± 

3 

0.1915 

-1 

10 

15 

-0. 1738 

-0.  1562 

10 

1 

± 

3 

0. 1887 

-1 

1  0 
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Table  07 


Zeros  of  ?irch/DS 


Speed 

Gain 

Boot 

1 

-4 

-0.4184  10 

-2 

-0.4370  10 

2 

-3 

-0.1674  10 

-2 

-0.8739  10 

3 

-3 

-0.3766  10 

-1 

-0.1311  10 

4 

-3 

-0.6675  10 

-1 

-0.1745  10 

5 

-2 

-0.1044  10 

-1 

-0.2182  10 

6 

-2 

-0.15C3  10 

-1 

-0.2619  10 

7 

-2 

-0.2047  10 

-1 

-0.3056  10 

8 

-2 

-0.2674  10 

-1 

-0.3493  10 

9 

-2 

-0.3385  10 

-1 

-0.3929  10 

10 

-2 

-0.4179  10 

-1 

-0.4367  10 

1  1 

-2 

-0.5057  10 

-1 

-0.4804  10 

12 

-2 

-0.6013  10 

-1 

-0.5238  10 

13 

-2 

-0.7058  10 

-1 

-0.5675  10 

14 

-2 

-0.8187  10 

-1 

-0.6 112  10 

15 

-2 

-C.9400  10 

-1 

-0.6549  10 
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Table  C8  -  Zeros  of  Pitch/DB 


Speed 

Gain 

Hoot 

1 

0.9111 

-5 

10 

0.9434 

-3 

10 

2 

0. 3644 

-4 

10 

0. 1887 

-2 

10 

3 

0.8200 

-4 

10 

0.2831 

-2 

10 

4 

0.  1453 

-3 

10 

0.3768 

-2 

10 

5 

0.2272 

-3 

10 

0.4713 

-2 

1  0 

6 

0.3274 

-3 

10 

0.5654 

-2 

10 

7 

0.4457 

-3 

10 

0.6599 

-2 

10 

8 

0.5822 

-3 

10 

0.7542 

-2 

10 

9 

0.7370 

-3 

10 

0 .8484 

-2 

1  0 

10 

0.  9  100 

-3 

10 

0.9428 

-2 

10 

1  1 

0. 1  10  1 

-2 

10 

0.1037 

-1 

10 

12 

0.1309 

-2 

10 

0.1131 

-1 

10 

13 

0. 1537 

-2 

10 

0.1225 

-1 

10 

14 

0.1783 

-2 

10 

0.1319 

-1 

10 

15 

0.2047 

-2 

10 

0.1414 

-1 

10 
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Table  09  -  Zeros  of  Depth/DS 


Speed 

Gain 

Root 

Root 

1 

-2 

-0.1904  10 

-2  -1 
-0.7732  10  ±  j  0.5083  10 

2 

-2 

-0.7618  10 

- 1  -  1 
-0. 1546  10  ±  j  0.4378  10 

3 

-1 

-0.1714  10 

-  1  -  1 
-0.2320  10  ±  j  0.2841  10 

4 

-1 

-0.3036  10 

-  1 

0.5803  10 

-2 

-0.3740  10 

5 

-1 

-0.4749  10 

-1 

0.9089  10 

-  1 

-0. 1366  10 

6 

-1 

-0.6643  10 

0.1183 

-1 

-0.2557  10 

7 

-1 

-0.9316  10 

0. 1438 

-  1 

-0.3567  10 

8 

-0. 1217 

0. 1685 

-  1 

-0.4489  10 

S 

■  -0.1541 

0. 1926 

- 1 

-0.5356  10 

10 

-C. 1902 

0.2165 

-1 

-0.6190  10 

1 1 

-0.2302 

0.2400 

-1 

-0.7001  10 

12 

-C.2737 

0.2633 

-1 

-0.7790  10 

13 

-C.3213 

0.2866 

-  1 

-0.8572  10 

14 

-0.3726 

0.3097 

-1 

-0.9343  10 

15 

-094278 

0.3329 

-0.1011 
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Table  10  -  Zeros  of  Depth/DB 


Speed 

Gain 

Root 

Root 

1 

-2 

-0.1106  10 

-1  -  1 
-0.  1030  10  ±  j  0.4792  10 

2 

-2 

-0.U422  10 

-1  - 1 
-0.2059  10  ±  j  0.4404  10 

3 

-2 

-0.995C  10 

-  1  -  1 
-0.3089  10  ±  j  0.3665  10 

-1 

-0.1764  10 

-  1  -  1 
-0.411210  ±  j  0.2269  10 

5 

-1 

-0.2757  10 

-1 

-0.7499  10 

-1 

-0.2784  10 

6 

-1 

-0.3972  10 

-0.  1049 

-1 

-0.1853  10 

7 

-1 

-0.54C6  10 

-0.  1304 

-1 

-0.  1359  10 

8 

-1 

-0.7C65  10 

-0.  1544 

-1 

-0.1020  10 

9 

-1 

-0.8944  10 

-0.  1775 

-2 

-0.7620  10 

10 

-C.  1  104 

-0.2003 

-2 

-0.5510  10 

1  1 

-0.  1336 

-0.2227 

-2 

-0.3718  10 

12 

-C.  1589 

-0.2447 

-2 

-0.2156  10 

13 

-0.  1865 

-0.2667 

-3 

-0.7430  10 

m 

-0.2163 

-0.2886 

-3 

0.5368  10 

15 

-0. 2484 

-0.3088 

-3 

0.1734  10 
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B.  STEADY-SIATE  DECOUPLING 


The  derived  lodel  represents  a  linear,  multivariable 
system.  In  order  to  allow  the  application  of  existing 
single-loop  techniques,  such  as  Nyquist-Bode-Nichol's 
methods  and  Root-Locus  design,  the  equations  must  be 
decoupled.  Decoupling  requires  the  system  to  be 
characterized  by  a  non- sing ular,  diagonal  transfer  function 
matrix.  >lhen  linear  state  variable  feedback  is  applied,  one 
necessary  condition  according  to  Wolovich  (Ref.  12)  is  that 
a  system  characterized  by  a  (n  *  m)  transfer  function  matrix 
may  not  have  poles  at  the  origin.  The  developed  linear 
system  contains  two  transfer  functions  with  such  poles. 


Ref.  4  develops  another  approach  which  decouples  the 
steady-state  only.  This  approach  uses  classical  cascade 
ccmpensaticn  and  allows  the  precluded  poles  at  the  origin. 
Actually,  these  poles  are  used  for  the  decoupling  of  linear, 
multivaria tie  systems. 


A  steady-state  decoupled  system  is  a  system  in  which 
changes  in  each  input  are  reflected  in  a  corresponding 
output  and  only  that  output,  when  steady-state  is  reached. 
Thus  mutual  interactions  are  allowed  during  the  transient 
period,  but  only  during  this  period. 


Using  cascade  compensation  and  a  diagonal  compensator 
matrix 


G  = 
c 


cl  1 


0 


0 


g 


c22 
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the  model  coEf iguration  becomes: 


where 


G 

P 


g  g 

?n  pi2 

g  g 

p21  p22 


By  theorem  2.2  of  Ref.  4  the  system  is  steady-state 
decoupled  if  and  only  if 


and 


lim 

s-*o 


1/  s 


(k-1) 


♦  (I+G  G  )  /det(I+G  G  ) 
p  c  1 2  pc 


=  0 


(5) 


2. 


(k-1) 

lim  1/  s  =»  (I  +  G  G  )  /det(I  +  G  G  ) 

s-*o  p  c  21  pc 


0 


(10) 
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since 


G  G 
P  c 


9  g  9  5 

pll  c11  p12  c22 


9  9  9  9 

p21  c11  p22  c22 


we  have 


det  (I+G  G  ) 
p  c 


1+9  9+99  + 

pll  cl  1  p22  c22 

9  9  9  9  “ 

p  1 1  p22  c 1 1  c22 

9  9  9  9 

pl2  p21  cll  c22 


and  cofactors 


(I  +  G  G  ) 
p  c  12 

(I+G  G  ) 

p  c  2 1 


9  9 

p21  cll 

9  9 

p  12  c22 


(11) 


(12) 

(13) 


K  and  k  are  defined  according  to 
1  2 


(input^ (t) )  =  1/  s 


where  the  inputs  are  constants,  steps, 
arbitrary  amplitude. 


ramps  etc, 


with 


By  substituting  (11)/  (12),  and  (13)  into  (9)  and 

(10)  ,and  assuming  a  step  input  for  depth  and  constant  pitch, 
we  have  the  two  necessary  and  sufficient  conditions  for 
steady-state  decoupling  : 


1 .  lim 

9  9 

/ 

D  =  0 

(14) 

s-»o 

p21  cll 

2 .  lim 

s*g  g 

/ 

o 

II 

o 

(15) 

S-^O 

12  22 

no 


4 

-  •  «  •  • 


i 


ifO  L»4 


.  -  .nifi. 


i';«| 


«< 


#  t  •« 
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where 


D  =  det  (I  +  G  G  ) 
F  c 


All  g  are  known.  Thas  the  design  for  steady-state 
F 

decoupling  is  the  determination  of  g  and  q  such  that 

cll  c22 

(14)  and  (15)  are  satisfied. 

Although  the  locations  of  the  poles  and  zercs  are 
dependent  on  the  speed  cheir  number  is  independent. 
Therefore  the  poles  will  oe  called 

p  DQ”  1^2^3... 

i  jm 

and  the  zeros 


z  0=  1,2,3... 

i  jm 

th 

where  ij  denotes  the  element  of  the  matrix  and  m  the  m 


pole  and  zerc  of  the  element. 


g 


g 


g 


g 


pii 

p12 

p21 

p22 


k  (s+z  )  (s+z  ) 
11  111  112 


k  (s  +  z  )  (s  +  z  ) 

12  121  122 


'2, 


k  (s  +  z  ) 

221  ^ 


/ 

/  sr 


where 


r= 

when  these  transfer  functions  are  substituted  into 
equations  (14)  and  (15),  both  limits  go  to  zero  for  all 
speeds . 

Thus  conditions  (1)  and  (2)  are  satisfied  fox  the  specified 
inputs. 

•  S§§i9Il  9I  SY  stem 

The  previous  section  resulted  in  the  decoupling  of 
the  steady-state  of  the  system.  However,  one  must  realize, 
that : 

1.  the  result  does  not  guarantee  stability. 

2.  steady-state  decoupling  is  only  meaningful, when 
the  closed-loop  system  is  stable. 

Thus  It  is  necessary  to  consider  the  stability  of  the 
clcsed-locp  system  after  the  steady-state  decoupling  has 
been  achieved. 

clef. 4  investigates  a  connection  between  single-loop 
and  multivariable  systems  by  properly  factorizing  the 
closed-lccp  characteristic  equation.  It  allows  the  design  of 
multivariable  systems  by  using  any  suitable  single-loop 
method.  The  primary  method  chosen  will  be  the  roct-locus 
design. 
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2.  Stability  criterion 

Ref.  13  and  14  prove  that  the  stability  of  this 
multivariable  system  is  determined  by  the  zeros  of  H^(s) 

and  N  (s) ,  where 
2 

N  (s)/D  (s)  =  det(I  +  G  {s)*G  (s))  (16) 

11  pc 

and 

N  (s)  =  A  c(s)  *Ap(s)/D  (s)  (17) 

2  1 

N^(s)/D  (s)  is  in  irreducible  fora,  i.e.  all  coamon  factors 

are  cancelled.  A c  (s)  represents  the  characteristic 

polynomial  cf  the  transfer  function  matrix  G  (s)  and  Ap(s) 

c 

represents  the  polynomial  of  the  transfer  function  matrix 

G  (s) .  The  characteristic  polynomial  of  the  transfer 
P 

function  is  defined  as  the  least  common  denominator  of  all 
miners. 


However  Ref. 4  devises  a  method  by  which  it  is 
necessary  to  check  equation  (16)  alone  for  stability,  if: 

1.  cancellations  are  selected  systematically  by 
using  equaticn  (4.5)  Ref.  4. 

2.  poles  of  G  are  carefully  selected,  which  can 

c 

be  taken  care  of  in  the  process  of  the  design. 
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Thus 


det  (I  +  G  G  )  =  0  (18) 

P  c 

represents  the  characteristic  eguarion  of  the  multivariable 
system. 

Fcr  the  system  under  consideration  equation  (18) 

becomes: 


det  (I  +  G  G  ) 
F  c 


1  +g  9 
pi  1  c  1 1 


g  g 

p22  c22 


det (G  ) 
P 


g  g 

cl  1  c22 


which  can  be  factored 


det  (I  +  G  G  )  =  (1 

P  c 


*  5  ,,9 

p  1  1  c  1  1 


+  G  g  ) 
eq  c22 


where 


G 

eq 


g 


(1  +  detG  *g  /g 

p  cl  1  p2  2 


)/(i+g 

pi  1 


(19) 


Ey  cancellation  of  the  common  roots  the  equation 
determining  the  stability  of  the  system  becomes 


1  +  G  g  =  0 
eq  c2 2 


(20) 


Equation  (20)  still  contains  both  g  and  g  as 

c 1 1  c22 

unlcnown  functions.  It  is  possible  for  example  to  choose  g 

c  1 1 

arbitrarily.  In  this  case  the  only  unknown  left  to  be 
designed  is  g  ,  which  means  the  design  is  reduced  to  a 


sicgle-lccp  case.  But  if  g  has  not  bean  chosen  properly, 

c  1  1 

G  may  become  very  unstable.  This  can  make  the  design  of 
eg 
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very  difficult.  Thus  it  is  preferable 


to  have 


c22 

guidelines  in  choosing 
pole-zero  pattern  for  G 


eg 

q  much  easier. 
c22 


g  in  order  to  get  a  reasonable 

cl  1 

which  will  make  the  design  of 


3-  Assign  Guidelines 

The  roots  of  equation  (19)  are  determined  by  the 
following  expressions: 


1-  =  g  +  det  (G  ) g 

z  p22  pell 


2 .  c-  =  1  +  g  g 

1  pll  cll 


where 


G  =  G  /G  (21) 

eg  1  2 


Then 


d€t(i  +  GG)  =  G  +  Gg 

p  c  1  2  c22 


As  the  design  is  concerned  with  the  poles  and  zeros 

of  G  they  must  be  identified.  They  can  be  found  by 
eg 

considering  the  functions  G  and  G  .  From  equation  (21)  it 


can  be  seen,  that  the  zeros  of  G  are  determined  bv: 

eg 


Z  =Z  +P  -(P  nP) 

Geg  G1  G2  G2  G1 


(22) 
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F  =Z  +P  -(P^P) 

Geq  G2  G1  G2  G1 


(23) 


Poles  and  zeros  of  G  and  G  are  dependent  on  the  function 

1  2 

q  acccrdinq  to: 
cll 


1  *  c€t(G  )*g  /  g  =  0 

p  cll  p  22 


(24) 


1  +  q  g  =0 
pi  1  cll 


(25) 


Therefore  the  design  will  proceed  as  follows: 


1.  Design  g  to  achieve  a  reasonable  pole-zero  pattern 
c  1 1 

for  G 

eg 


a.)  Prepare  root-locus  studies  for  equation  (24)  and 

(25),  where  g  =  f (K  ,s)  .  K  is  a  real  parameter 

cll  1  1 

positive  or  negative.  The  pcles  and  zeros  of  g 


are 


c  1  1 


determined  by 


aa. )  number  of  free  integrators  needed  to  decouple 
the  system. 

ab. )  the  needs  of  compensation,  which  may  arise, 
due  to  desired  root  locations. 

b.)  Frci  root-locus  studies  choose  the  parameter  K^, 

which  fixes  g  .  This  will  be  an  iterative  process  as 
cl  1 

both  rcct-lccus  studies  in  combination  with  equations 

(22)  and  (23)  determine  the  final  pole-zero  pattern  of 

G  « 

eq 


1 16 


1  L  • 


I 

« 


I 


2.  Design  cf  g  in  order  to  meet  the  system 

c22 

specification.  As  G  is  icnown  equation  (20),  i.e. 

eg 


1  +  G  g  =0 
eg  c22 

can  be  studied  via  root-locus,  where  g  =  f (K  , s) .  The 

c2  2  2 

criteria  which  determine  f (K  , s)  are  the  same  as  mentioned 

2 

above  for  f  (K  ^  ,s)  . 

3.  Wirh  the  proper  choice  of  K 

2 

a. )  system  stability  is  guaranteed. 

b. )  the  dynamics  of  the  system  are  fixed. 

c. )  the  diagonal  compensator  is  determined. 

4.  The  design  is  completed.  Mow  simulation  is  required  to 

check  whether  the  design  specifications  have  been  met.  If 

the  specifications  have  not  been  met,  a  redesign  of 

compensator  g  or  even  g  may  be  necessary.  But  every 
c22  cl  1 

trial  provides  some  insight  and  helps  with  the  next  trial. 


^  design 


Following  the  design  guide  lines  the  first  step  is  a 

loot-locus  study  for  the  equations  (24)  and  (25).  As  shown 

in  II. B  no  free  integrators  for  g  and  g  are  required 

c11  c22 


for  the  specified  inputs. 


The  speed  range  of  interest  extends  from  3  kn  to  15 


•  i  1  3 


1 


I  t 
t  .  . 


•  AT  ji 

•  n  *3*41 

»« 

^  .til 


I 


i  >Oi 


4 


kn.  As  the  root  locations  of  the  system  change  with  varying 
speed,  it  is  not  certain  whether  one  compensator  set  will 
cover  the  whcle  speed  range.  It  seems  reasonable  to  design 
for  the  medium  speed  of  9  kn.  Afterwards  it  can  be  tested 
whether  this  compensator  set  is  adequate  for  the  whole  speed 
range  and  whether  the  variations  due  to  speed  can  be  taken 
care  of  by  either  gain  variations  or  other  means.  It  may 
turn  out  that  different  compensators  are  required  tc  meet 
the  design  specifications. 

a.)  Design  of  g 

c  1 1 

The  rcct-lcci  for  1+g  g  =0  where  g  =K  are 

c n  pii  ell  1 

shewn  in  Fig.  70  K  =  negative  and  Fig.  71  K  =  oositive. 

1  21  ‘ 

=  positive  is  clearly  undesirable,  because  there  will  be 
one  branch  that  extends  along  the  positive  real  axis  tc  + 
and  will  produce  a  pole  on  the  positive  real  axis  for  G 

eg 

The  root-loci  for  1  +  det  (G  )  *g  /g  =0  are 

p  c  1  1  p22 

shewn  in  Fig.  72  K  =  negative  and  rig.  73  K  =  positive. 
For  the  same  reason  as  above  has  to  be  negative. 

From  the  root-loci  for  K  =  negative  it  is  obvious 

1 

that  further  compensation  is  needed  in  order  to  achieve 

reasonable  root  locations.  The  following  compensator  was 
chesen 


2  2 

c  1 1  1 


118 


.1  UNITS  PER  INCH 

Figure  70  -  ROOT  LOCOS  1  +  GPn*GC11  =  0 
SPEED  =  9  KN,  K1  =  NEG . 
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.1  UNITS  PER  INCH 

Figure  71  -  ROOT  LOCOS  1+G?11*GC11  =  0 
SPEED  =  9  KN,  K1  =  POS. 
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.1  UNITS  PER  INCH 


Figure  72  -  ROOT  LOCOS  1 > DET ( GP) /GP22*GC1  1  =  0 
SPEED  =  9  KN,  K1  =  NEG. 
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.1  UNITS  PER  INCH 


Figure  73  -  ROOT  LOCUS  1+DET (GP) /GP22*GC1 1  =  0 
SPEED  =  9  KN,  K1  =  POS. 
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This  compensation  resulted  in  Pig.  74  for 

1  +  c  g  =0  and  Fig.  75  for  1+det  (G  )g  g  =0.  The 
'c11  p11  p  c11  p11 

rocts  of  Fig.  74  represent  part  of  the  poles  of  G  -By 

eg 

means  of  the  compensation  all  roots  were  constrained  tc  the 

stable  region.  The  roots  of  Fig.  75  will  form  the  zeros  of 

G  .  Although  this  compensation  did  not  pull  the  unstable 
eg 

branches  considerably  to  the  left  half  plane,  it  does  not 

mean  that  the  system  will  not  be  well  behaved,  because  the 

system  behavior  will  in  the  end  be  determined  by  the  poles 

of  G  and  g 

eg  c22 

As  Sef.  4  in  eguations  (6-33)  points  out,  the  poles 

and  zeros  cf  G  are  determined  as  follows: 
eg 


Poles:  P 

=  Z 

+  P  - 

(?  ^ 

P 

Geg 

G1 

G2 

G2 

G1 

Zeros;  Z 

=  Z 

+  P  - 

(P  n 

P 

Geg 

G2 

G  1 

G2 

G  1 

P  =  P  +  p 

G1  gpll  gcll 


) 

) 


Z  =  frcots  from  foot-locus  1+g  g  =0) 
G1  p11  c11 


P 

G2 


P 

gp22 


+  P  (detG  )  g  -  (P 

p  ell  gp22 


/^  P(detG  )  g  ) 
p  c  1 1 


B  ^  P  =  P  for  this  system. 

G2  G1  G1 
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.  1  UNITS  PER  INCH 


Figure  74  -  ROOT  LOCUS  1  +  GPl1*GCn  = 
SPEED  =  9  KN,  K1  =  NEG. ,  COMPENSATED 
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.1  UNITS  PER  INCH 


Figure  75  - 
SPEED  = 


ROOT  LOCOS  1+DET (GP) /GP22*GC1 1 
9  XN,  K1  =  NEG.,  COMPENSATED 


0 


125 


Thus  the 


fcles  of  G 


are : 


P 

Geg 

and  the 

Z 

Geg 

system  t 
zeros  h 
studies 

behave  w 
Fig.  76 
They  are 

b.)  Desi 

plane  d 
only  to 
right  g 

determin 

G  ,  thu 
eg 


eg 


P  -{0.0}+  (roots  from  root-locus  1+g  g  =0} 
gpll  pllcll 


zeros  of  G  are: 
eg 


(roots  from  root-locus  1+detG  g  /g  =0}  . 

p  c 1 1  p22 


G  has  four  known  poles,  namely  the  ooles 
eg 


ransfer  fu 

notion  g 

P 

1  1 ' 

The 

remaining 

pole 

ave  to  be 

deter  mined 

from 

the 

compensated 

roct 

Fig.  74  a 

nd  Fig.  75 

cf  the 
s  and 
-locus 


Ic  get  an  idea  how  the  pole-zero  pattern  of  G  will 

eg 

ith  varying  both  root-loci  are  superimposed  in 

The  additional  poles  of  g  -  {0.0}  are  added. 

pll 

independent  of  K  and  therefore  fixed. 

1 


gn  cf  g 

c22 
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.1  UNITS  PER  INCH 


Figure  76  -  SUPERPOSITION  OF  FIG.  74  AND  FIG.  75 
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ooen-loop  system  G  .  The  closed-loop  poles  are  now 

eg 

determined  by  the  root-locus  study  for  1+G  g  =0.  Where 

eg  c22 

q  =  K  .  The  gain-compensated  studies  for  K  =1.48  and  K 
c22  2  1  1 

=  75.0  are  shown  in  Fig.  77  and  Fig.  78  .  Fig.  77  has 

complex  roots  clcse  to  the  imaginary  axis  which  go  to 

asymptotically,  which  is  not  desirable  as  the  damping  gets 

too  small.  Fig.  78  shows  very  similar  behavior.  Therefore 

it  is  desired  to  move  the  asymptotes  to  the  left,  which  is 

achieved  by 


g  =K  (s+0.5)/(s+  10.0)  . 
c22  2 

The  results  of  this  compensation  are  shown  in  Fig.  79 
Fig.  80.  for  K  =  1.48  and  K  =  75.0  respectively.  As 


compensated  root- locus  for  K 
stable  for  all  values  of  K  *K 


1  .48 
The 


indicates , 
root- Iccus 


G  will 
eg 

for  .K 


2  eg  1 

75.0  has  zeros  in  the  right  half  plane  and  will  limit 

value  of  K  *K  in  order  to  maintain  stability. 

2  eg 


and 

the 

be 

the 


As  many  poles  are  clustered  around  the  origin,  it  is 
difficult  to  tell  which  roots  will  be  dominant.  A  time 
response  could  reveal  the  dominant  roots,  but  as  the  system 
has  to  be  simulated  in  order  to  check  the  compensation 
anyhow,  because  both  outputs  have  to  satisfy  the  design 
specifications,  it  seems  more  reasonable  to  simulate  the 
total  system.  From  the  simulation  results  it  will  be  judged 
whether  the  compensation  was  sufficient  or  if  further 
compensaticD  is  necessary. 
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.  1  UNITS  PER  INCH 


Figure  77  -  ROOT  LOCUS  1+G(ZQ)*GC22  =  0 
SPEED  =  9  KN,  K1  =  1.48 
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Figure  78  -  ROOT  LOCUS  1 +G  (EQ)  *GC22  =  0 
SPEED  =  9  KN,  K1  =  75. 
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.1  UNITS  PER  INCH 


Figure  79  -  ROOT  LOCUS  1+G(EQ)^GC22  =  0 
SPEED  =  9  KN,  K1  =  1.48,  COMPENSATED 
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Figure  80  -  ROOT  LOCUS  1+G(EQ)*GC22  =  0 
SPEED  =  9  KN,  K1  =  75.,  COMPENSATED 
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5. 


2^  compensated  linear  model 


a.  Determination  of  gain  constants 

As  mentioned  above  the  roots  of  G  will  always  be  stable 

eg 

if  a  value  of  K  =  1.48  is  chosen.  Thus  this  value  was 

1 

picked  for  the  initial  simulation  trials.  The  root-locus  for 

G  for  this  choice  of  K  is  shown  in  Fig.  79  .  The  next 
eg  1 

step  is  to  find  the  gain  constant  K^.  The  gains  shown  in 

Fig.  79  represent  K  ♦K  .  In  order  to  determine  K  ,  K 

2  eg  2  eg 

has  to  be  found. 

Because 


G 

eg 


g  +  de t  (G 

p22 


)  *9 
P 


c22 


/1+g 


pll 


g 


cl  1 


-2  7  9 

=  (-.  3305  1  0  )  (s  +...)/{s  +...) 


-2 

K  =  -.3385  10 
eg 


Therefor  e 


K  =  (gain  from  root-locus) /K 
2  eg 

The  only  roots  which  change  significantly  with  increasing 
gain  are  those  which  go  to  .  A  reasonable  value  for  K 


60,  because  ; 


seems  to  be  K  = 
2 


-  K  is  not  too  large  (to  avoid  excessive  rudder 


deflections) 


-  damping  is  greatest  in  this  region. 


Six  siiulaticn  runs  were  performed  using  the 
as  pointed  out  above  and  variatons  of  them, 
depth  change  is  10  ft.,  the  ordered  pitch  is  0 
The  results  are  shewn  in  Table  11. 


chosen  values 

The  ordered 
o 
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Table  11  -  Parameter  study  for  gains 
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Figure  81  -  RON  1 
10  FT  DEPTH  CHANGE,  9  KN 
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Figure  82  -  RON  2 
10  FT  DEPTH  CHANGE,  9  KN 
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Figure  83  -  RON  3 
10  FT  DEPTH  CHANGE,  9  KN 
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Figure  84  -  RON  6 
10  FT  DEPTH  CHANGE,  9  KN 
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Figure  85  -  RON  6 
10  FT  DEPTH  CHANGE,  9  KN 
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Figure  86  -  RUN  6 
FT  DEPTH  CHANGE,  9  KN 
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All  these  test  runs  showed  that  the  compensated  system  is 
stable  as  expected  and  that  the  steady-states  are  reached 
independently  without  error.  The  plane  angles  increase  with 
increasing  gains.  The  pitch  behavior  has  basically  the  same 
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Run  1  : 

K  =  1 .  K  =15. 

1  2 

The  ordered  depth  was  reached  with  very  small  overshoot,  but 
the  system  showed  a  big  undershoot  of  10.8-;?  after  the  first 
peak.  This  is  due  to  the  fact  that  at  the  time  the  ordered 
depth  is  reached,  the  pitch  is  greatest  in  the  positive 
direction.  This  is  not  desirable. 
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Run  2  : 


K  = 

1 

The  increase  in 
plane  deflection. 


.5  K  =  15. 

2 

K  leads  to  an  increase  in  the  fairwater 

1 

The  depth  change  is  accelerated  and  the 


pitch  deviations  increase,  which  in  turn  leads  rc  greater 

stern  plane  deflections.  The  amplitudes  of  the  oscillations 

increase,  the  overshoot  reaches  11.4%  and  the  undershoot 

11.4%.  As  the  fairwater  planes  have  reached  about  the 

allowed  deflection  limits  K  will  not  be  increased  anv 

1 


further.  The  action  taken  will  concern  K  onlv. 

2 


Run  3  ; 


K  =  1 .5  K  =  40. 

1  2 

The  increase  in  resulted  in  larger  stern  plane 

deflections,  which  damped  the  pitch  more  than  before.  Due  to 

the  smaller  pitch  changes  which  at  the  same  time  support  the 
depth  response,  the  deviations  are  decreased  significantly 
and  die  out  rapidly.  The  trade-off  which  has  to  be 
considered  is  that  increasing  the  damping  will  result  in 
longer  tines  until  the  ordered  depth  is  reached. 

Run  4  and  5  : 

are  not  considered  as  mentioned  above. 
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Bud  6  : 


K  =1.5  K  =60. 

1  2 

This  further  increase  in  K  confirmed  the  previously 

2 

observed  tendency.  The  overshoot  is  decreased  to  3.6  %  and 

there  is  essentially  no  undershoot.  Both  run  3  and  6  satisfy 
the  design  specifications.  But  run  6  shows  clearly  the  more 
desirable  behavior  as  the  pitch  is  smaller,  the  overshoot  is 
less  and  the  depth  accuracy  is  more  important  than  the  small 
sacrifice  in  time,  as  all  these  depth  maneuvers  represent 
standard  depth  changes  and  no  emergency  dives. 


b.  Check  for  valid  speed  range 


So  far  the  compensator  is  designed  for  9  kn.  The 
sensitivity  cf  the  compensated  system  to  speed  is  important 
for  two  reasons.  First,  the  constant  speed  assumption  for 
the  linear  model  is  not  quite  accurate,  because  on  our 
submarines  the  speed  is  ordered  in  fora  of  rpm.  The  actual 
speed  varies  with  rudder  deflections  and  pitch  angle.  These 
influences  result  in  small  speed  changes.  This  phenomenon 
was  observed  in  the  comparative  test  runs  with  the  linear 
and  non-linear  models.  Second,  it  is  of  interest  to  find  the 
speed  range  where  this  compensator  still  yields  a 
satisfactory  behavior  which  is  within  the  design 
specifications. 


designed 
kn  to 
Table  12 


Therefore  simulation  runs  with  the  compensator 
for  9  kn  were  performed  for  the  speed  range  from  3 
15  kn  at  2  kn-intervals .  The  results  are  shown  in 
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Table  12  -  10  ft  depth  change 
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Figure  89  -  RON  3 
100  FT  DEPTH  CHANGE,  7  KN 
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Figure  91  -  RON  6 
100  FT  DEPTH  CHANGE,  13  KN 
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Figure  92  -  RUN  7 
100  FT  DEPTH  CHANGE,  15  KN 
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As  can  be  seen  from  Table  06  the  poles  of  the  characteristic 
equation  of  the  open  loop  system  transfer  functions  move 
closer  tc  the  origin  of  the  s-plane  with  decreasing  speeed, 
which  should  lead  to  a  slower  response.  At  the  same  time  the 
damping- ratic  decreases  with  decreasing  speed.  All  test  runs 
reflect  this  observed  tendency. 


*  Run  4-7  ; 

The  time,  at  which  the  ordered  depth  was  reached  for  the 
first  tine,  is  acceptable.  The  overshoot  is  within  the 
design  specifications.  The  plane  angles  never  exceeded  the 
limits  and  the  pitch  is  small.  Small  gain  variations  led  to 
a  shift  in  the  trade-off  between  response  time  and 
overshoot.  The  improvements  which  resulted  from  tests  were 
not  significant  enough  to  justify  the  complexity  of  gain 
switching.  Thus  the  designed  compensator  will  cover  the 
entire  speed  range  from  9  to  15  kn  very  satisfactorily. 


*  Run  3  : 


The  7  kn  depth 
slightly  with 
plane  angles 
represents  a 
acceptable . 


response  exceeds  the  design  specifications 
respect  to  time  and  overshoot,  while  pitch  and 
are  still  within  the  limits.  This  speed 
borderline  case  where  the  compensator  might  be 


Run  2  and  1  : 


For  these  two  speeds 
respectively  exceed 
required  are  clearly 


the  overshoot  of  11.58  %  and 
the  specifications  by  far  and 
not  desired. 


18.77  % 

the  times 


It  was  possible  to  improve  the  response  time 
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with  gain  variations,  but  it  was  impossible  to  reduce  the 
overshoot  to  an  acceptable  value.  This  indicated  that  for 
slew  speeds  gain  switching  can  not  be  a  sufficient  measure 
to  adjust  this  compensator.  It  seems  that  a  redesign  of  the 
compensatcr  is  necessary. 


3  kn  and  5  kn  are  eliminated  from  further 


considerations  as  these  speeds  are  not  too  important  for 
submarine  maneuvering. 


c.  Modification  for  major  depth  changes 


The  mechanical  limit  for  both  rudder  deflections  is  35  °. 
AS  outlined  in  the  design  specifications  these  limits  are 
reduced  to  smaller  values  for  higher  speeds.  Due  to  the 
linearity  of  the  system  major  depth  changes,  which  lead  to 
proportional  depth  errors,  will  cause  plane  deflections 
exceeding  these  limits.  Therefore  it  is  necessary  to  prevent 
the  planes  from  exceeding  these  limits.  There  are  two 
possible  actions  one  can  take. 

First,  limit  the  planes.  Test  runs  where  this 
scheme  was  implemented  led  to  unacceptable  plane  behavior. 

Second,  limiting  of  the  error  signals.  The 
resultant  block  diagram  is: 
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From 


the  design  specifications  it  was  determined  that  the 
error  limit  function  could  be  approximated  by  a  straight 
line ; 

lim  =  -1.6655*ac  +  59.6875 

I 

Where  UC  is  in  ft/s. 

This  function  and  the  resulting  maximum  angles  are  shown  in 
the  following  diagram 


Table  13  shows  the  simulation  runs  for  100  ft  depth 


change. 
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Table  13  -  100  ft  depth  change 
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Figure  93  -  RON  1 
100  FT  DEPTH  CHANGE,  7  KN 
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Figure  94  -  RON  2 
100  FT  DEPTH  CHANGE,  9 
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Figure  96  -  RUN  4 
100  FT  DEPTH  CHANGE,  13  KN 
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Figure  97  -  RON  5 
100  FT  DEPTH  CHANGE,  15  KN 
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Figure  98  -  RUN  5 
100  FT  DEPTH  CHANGE,  15  KN 
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After  an  initial  pitch  the  su 
depth  with  constant  plane 
depth  is  reached,  opposite 
eliminate  the  depth  rate. 

Although  it  takes 
depth  initially  the  overall 
with  increasing  speed. 


baarine  approaches  the  ordered 
deflections.  Shortly  before  the 
rudder  angles  are  used  to 

longer  to  reach  the  ordered 
time  to  steady-state  decreases 


This  completes  the  design  for  the  speed  range 
from  7  kn  to  15  kn  and  all  possible  depth  changes.  From  here 
on  the  behavior  of  the  compensated  non-linear  model  will  be 
investigated  under  special  conditions,  like  turns,  out  of 
trim  conditions, etc. 

The  final  block  diagram  and  the  DSL-Program  are  shown  in 


App.  D. 
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C.  TOTAL  DECOUPLING 

The  previous  design  led  to  a  steady-state  decoupled  system, 
where  both  outputs  were  influenced  by  a  single  input  during 
the  dynanic  period.  In  our  case  the  pitch  encountered  was 
small.  In  order  to  keep  the  submarine  at  even  keel  for  zero 
pitch  input  at  all  times  the  system  has  to  be  totally 
decoupled.  One  method  to  achieve  this  behaviour  for  a 
multivariable  system  will  be  analysed  and  applied  to  the 
linear  mcdel. 

The  system  configuration  chosen  for  this  design  is: 


The  closed  loop  system  can  be  represented  by; 


which  can  be  combined  to,  as  will  be  shown: 


U 
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T  is  the  desired  overall  transfer  function  which  determines 
the  desired  overall  respones. 

Therefore 

T  =  FG 


thus 


G  = 


-1 

F  T 


Frcm  the  block  diagrams: 

Y  =  TO 

Y  =  FX 

X  =  CE 

E  =  0-FY 


Thus 


Y  =  FC[U-FY] 


or 


-  1 

Y  =  [I  +  PCF ]  PCU 


and 


-  1 

T  =  [I  +  PCF]  PC  =  PG 
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Figure  100  -  BLOCK  DIAGRAM 
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As  can  b€  seen  from  the  block  diagram  the  requirements 
for  total  decoupling  can  be  stated  in  a  different  way. 

It  is  required  that 


eg  =  -c  g 

11  p21  21  p22 


and 


(26) 


eg  =  -c  g 
12  c11  22  p12 


(27) 


Single  loop  methods  can  be  applied  to  design  the 

compensators  c  and  c  to  guarantee  the  desired  response. 

1  1  22 

Once  c  and  c  are  determined  the  compensators  c  and  c 

11  22  12  21 

can  be  calculated  from  (26)  and  (27)  . 


-1 

As  Ref.  15  shows  C  is  realizable  if  G  and  [I  -  T]  are 
realizable.  A  rational  function  in  s  is  physically 
realizable  as  a  vcltage  transfer  function,  if  the  following 
conditions  are  met: 


the  coefficients  of  the  function  must  be  real, 
the  order  of  the  numerator  must  be  less  or  equal 
to  the  order  of  the  denominator. 

there  are  no  poles  in  the  right  half  plane  and  poles 
on  the  imaginary  axis  must  be  simple. 


1 .  Cesign 

It  is  desired  that  the  transient  response  has  no  overshoot 
and  has  an  acceptable  rise  time.  The  following  transfer 
function  is  chosen  for  both  outputs 
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2 


t  =  t 
1  1  22 


s  +«<. 


t  =  t  =0,  because  the  system  is  to  be  decoucled, 
12  21 


Ecllcwing  the  previously  outlined  design  steps. 


-1 

P  *T  has  to  be  determined, 


The  plant  matrix  has  the  following  elements 


=  k,, (zeros  gp,,)/ss- 


g  =  k  (zeros  g  )/sX" 

p12  12  p12  ® 


g  =  k  (zeros  g  )/v- 

p21  21  ^n7^''0 


p21 


g  =  k  (zeros  g  ) /^ 
p22  22  p22  ^ 


Thus 


-1 


(zeros  g  -k  (zeros  g  ^)/rs 

22  p22  0  12  p12  ® 

•k  (zeros  g  ) /vr  k  (zeros  g  )/v-s 

21  p21  0  11  ?1 1  ^ J 


g  g  -  g  g 

p11  p22  p12  p2l 


It  can  be  shewn  that 


g  g  ■  g  g  =  Vsr 

p11  p22  p12  p21 
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The  refore 


-1 


F 


k  (zeros  g  )s  -k  (zeros  g  ) 
22  p22  12  pl2 


(zeros  g  )  s 

p21 


k  (zeros  g  ) 
11  pi  1 


P 


2  1 


D 

‘  22 


and 


-1  2  2-1 
G  =  F  T  =  [  ol.  {s  *  ck)  ]p 


The  next  step  is  to  find 


[I-T  j 


-1 


l-r 


22 


1-t 


1  1 


(1-t  )  (1-t  ) 

1  1  22 


2 

(s+c<)  /s(s+2c^)  0 


2 

0  (s+d)  /s  (s+2d) 


/ 
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Finally  the  ccmpensator  matrix  C  can  be  determined; 


-1 

C  =  G[I-T] 


c  c 

11  12 


c  c 

21  22 


where 


c  =  oC.  k  (zeros  g  )/(s+2ol) 

11  2  2  p  22 


2 

c  =  -<^  k  (zeros  g  )/s(s  +  2o0 

12  12  p12 


c  =  -dk  (zeros  g  )/(s+2d) 

21  21  p21 


2 

c  =  k  (zeros  g  )/s(s  +  2o() 

22  11  p11 


The  elements  of  the  compensator  for  this  specific 
design  consist  of  the  zeros  of  the  plant  transfer  function 
and  poles  which  depend  on  the  desired  response.  As  the  zeros 
of  the  plant  and  probably  the  time  constants  of  the  desired 
response  depend  on  the  speed  of  the  submarine,  it  can  be 
expected  that  the  compensator  is  highly  dependent  on  speed. 

The  realizability  of  C  is  guaranteed,  if  G  and 

-  1 

[I  -  T]  are  realizable.  It  can  easily  be  seen  that  for  all 
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- 1 

elements  of  G  and  [I  -  T]  all  three  above  mentioned 
conditions  are  satisfied. 


2 .  Simulation 


In  order  to  check  the  design  the  system  was  simulated  for  a 
speed  of  S  kn.  The  choice  of  the  overall  transfer  function 


2  2 

t  =  t  =  .  2  /  (s+  .  2) 
11  22 


t  =  t  =  0 
12  21 


should  lead  to  an  overdamped  response,  which  reaches  the 
ordered  depth  in  about  60  sec.  From  Tables  06  -  10  and  the 
defined  overall  T  it  follows  : 


-1 

c^^  =  -.  31  945  10  (25.U52S+1)/ (2.5S+1) 


c^^  =  -. 38  194  (-5.  192  1S+1)  (18.671s>1) /s(2.5s+1) 


-2 

c^^  =  .  150  17  10  (-  1  17. 87S+1)  /  (2 .5s+1) 


-  1 

c^  =  -.  29056  10  (5.6338S+1)  (  131 .233+1)/s (2. 5s+ 1) 


*  Unconstrained  trials: 

In  the  first  simulation  the  submarine  had  to  change  the 
depth  by  10  ft  at  even  keel.  The  simulation  showed  the 
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expected  response,  i.e.  the  depth  was  reached  after  50  sec. 

-5 

without  overshoot.  The  pitch  was  of  the  order  of  ±10  .  This 

deviation  can  be  explained  by  round-off  and  truncation. 

In  the  design  the  physical  constraints  on  the  planes  were 
neglected.  As  it  turned  out  these  constraints  were  exceeded 
by  far.  Due  to  the  step  input  of  the  ordered  depth  the 
fairwater  plane  angle  increased  initially  to  117°  and 
exceeded  the  limit  of  30°  for  about  4.5  sec.  The  max.  stern 
plane  angle  was  about  25°.  Obviously  further  compensation 
is  necessary. 

Another  important  aspect  is  the  sensitivity  w.r.t. 
speed.  Therefore  the  submarine  was  run  with  the  same 
compensator  for  7  and  11  icn. 

For  7  kn  the  depth  response  was  oscillatory  with  an 

overshoot  of  1.6%.  The  pitch  oscillated  with  a  amplitude  of 

/ 

•  I  • 

For  11  kn  the  depth  response  was  oscillatory  with  an 
overshoot  of  IJ?.  The  pitch  oscillated  with  an  amplitude  of 
.09°.  As  expected  the  system  is  no  longer  decoupled,  because 
the  "correct"  compensator  depends  among  others  on  the  zeros 
of  the  system,  which  in  turn  vary  with  speed.  As  the  depth 
response  deviates  significantly  from  the  desired  response 
with  change  in  speed,  it  can  be  concluded  that  the 
compensator,  designed  for  a  specific  speed,  covers  only  a 
relatively  small  interval. 

*  Implementation  of  constraints  ; 

The  two  constraints  are  the  rate  limit  of  the  planes  and  the 
max.  angle  allowed  for  a  specific  speed.  The  rate  limit  is 
taken  care  of  by  the  actuator.  The  angle  can  be  kept  within 
its  limits  (ca.  30°  for  9  kn)  using  either  a  limit  on  the 
signal  applied  at  the  actuator  or  error  limitation.  For  a 
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Figure  101  -  TOTAL  DECOUPLING 
10  FT  DEPTH  CHANGE,  9  KN 
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Figure  102  -  TOTAL  DSCODPLING 
10  FT  DEPTH  CHANGE,  9  KN 
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Figure  103  -  TOTAL  DECOOPLING 
10  FT  DEPTH  CHANGE^  9  KN 
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depth  change,  generally  the  required  fairwater  plane 
deflections  are  much  larger  than  the  stern  plane 
deflections.  As  soon  as  only  one  plane  limiter  is  saturated, 
the  achieved  decoupling  is  disturbed  as  can  be  seen  from 
Fig.  100.  This  saturation  occurs  for  a  small  depth  change  of 
only  3  ft.  Obviously  this  scheme  is  not  acceptable. 

The  error  limitation  eliminates  these  problems.  When 
the  depth  error  is  limited  to  2.5  ft.  the  initial  fairwater 
plane  angle  stays  within  the  specifications.  On  the  other 
hand  the  "steady-state"  angle  for  large  depth  changes  turned 
out  to  be  toe  small  (4.5°) ,  which  leads  to  an  unacceptaole 
response  time.  To  achieve  both  an  acceptable  initial  and 
"steady-state"  angle,  the  limit  was  icept  at  2.5  ft  for  1.5 
sec  and  then  increased  to  a  final  value  of  12.5  ft  at  a 
constant  rate. 

The  limiting  function  used  is 


Lia  =  2.5  ft  t  <  1.5 

Llh  =  .2*  (TI«E-1 .5)  +2.5  ft  1.5  <  t  <  53.5 

LIM  =  12.5  ft  t  >  53.5 

The  final  block  diagram  and  the  DSL-Program  are  given  in 
App  .  D 

Fig.  101  to  106  show  a  10  ft  and  100  ft  depth 
change.  In  hoth  cases  the  depth  response  is  overdamped  and 
satisfies  the  time  requirements.  The  pitch  is  decoupled  and 
the  submarine  remains  at  even  keel. 


It  has  been  shown  that  this  design  approach  has  the 
fcllcwing  disadvantages: 


Speed  range  which  is  covered  by  individual 
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controller  is  very  small.  If  the  speed  deviates  too  much 
from  the  design  speed  the  system  is  no  longer  totally 
decoupled . 


Ihe  high  speed  dependence  requires  a  ccmplex 
controller  with  varying  pole-zero  combinations  dependent  on 
the  speed  variations  (  system  dynamics  ) . 

Additionally  the  "  Total  decoupling  scheme  ”  is 
disturbed  by  the  non-linear  terms  in  the  ncn-linear  model, 
which  are  not  contained  in  the  linear  approximation. 
Therefore  this  design  approach  will  not  be  pursued  further. 
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V.  SIMULATION  OF  THE  COMPENS^ED  NONrLIM45  ^OCSL 


A.  MODIFICATIONS  FOR  NON-LINEAR  MODEL 

Two  of  the  linearizing  assumptions  were,  that  the  forward 
speed  is  constant  and  that  the  submarine  is  always  in  trim. 

The  purpose  of  this  chapter  is  to  investigate  the 
deviations  from  these  assumptions  and  their  influence  cn  the 
behavior  cf  the  submarine  during  basic  maneuvers  like  depth 
changes,  turrs,  and  turns  during  depth  changes. 

Generally  trim  is  a  moment,  but  it  is  commonly  used  in  a 
different  context  ; 

-  "In  trim"  has  the  meaning  that  the  submarine  maintains 
depth  at  a  given  speed  with  the  desired  pitch  angle. 

-  "Out  of  trim"  conditions  are  heavy  or  light  forward, 
heavy  or  light  aft,  heavy  or  light  overall,  and  combinations 
of  these. 

-  '«To  trim"  means  to  change  the  contents  of  the  tanks  by 
shifting  water  between  tanks,  flooding  tanks  from  sea,  or 
pumping  water  from  tanks  ro  sea. 

The  arrangement  of  depth  control  or  auxiliary  and  trim  tanks 
in  the  assumed  submarine  are  such  that  any  change  in  the 
depth  control  tank  does  not  produce  a  moment,  i.e.  the  depth 
control  tank  is  at  the  center  of  gravity.  The  trim  tanks  are 
arranged  as  shown  in  Fig.  107. 


1 .  Trim 

Frcm  experience  it  is  known  that  a  submarine,  which  is  in 
neutral  conditions  for  a  certain  speed,  "  gains  buoyancy" 
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Figure  107  -  CONFIGURATION  OP  TANKS 
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WATER  FLOODED  INTO  AUXILIARY  TANK  IN  1000  LBS 


12 


SPEED  IN  KN 

Figure  108  -  CORHECTIONS  FOR  TRIM 
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WATER  TRIMMED  FORWARD  IN  1000  LBS 


and  pitches  up  with  increasing  speed  and  vice  versa  due  to 
dyramic  forces  acting  on  the  hull,  and  therefore  is  no 
longer  in  trim  as  defined  above. 

Generally  a  submarine  should  be  in  trim  at  any  time 
for  a  given  commmand  speed  to  avoid  initial  plane 
deflections.  This  guarantees  availability  of  the  full  range 
of  the  planes  for  maneuvering.  Although  it  is  nearly 
impossible  tc  get  the  exact  trim,  it  certainly  is  helpful  to 
knew  the  corrections  for  different  speeds, 

Assuning  neutral  conditions  for  0  kn  speed,  the 
corrections  for  trim  were  found  by  trial  and  error.  The 
values  shewn  in  Fig.  108  guarantee  an  approximately  neutral 
trim. 


The  results  confirmed  the  tendencies  observed  in 
practice . 


2‘  dfchanic^  constraints 

The  design  specifications  .  define  the  maximum  plane 
deflections  as  a  function  of  speed.  While  the  ordered 
revolutions  cf  the  shaft  remain  constant,  the  actual  forward 
speed  (  in  direction  of  the  x-axis  of  the  body-fixed 
coordinate  sjstem)  decreases  during  maneuvers  due  to  drag, 
which  is  produced  by  the  deflected  planes,  the  pitch,  the 
roll,  and  the  yaw. 

A  feedback  of  the  actual  forward  speed  U  instead  of 
the  ordered  speed  UC  modifies  the  error  limiter  to 


LIM  =  -1.6655*U  +  59.6375 


and  allows  the  use  of  the  maximum  plane  deflections  as  a 
function  cf  the  actual  speed. 
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TEST  EUNS 


For  the  test  runs  to  follow  the  designed  coa pensator , 
the  modified  limiter,  and  the  trim  values  for  the  ordered 
speed  according  to  Fig.  108  are  implemented  in  the 
non-linear  mcdel. 

The  program  and  final  block  diagram  are  given  in  App.  D. 


1  changes 

Table  14  lists  the  results  of  the  simulations  from 
7-15  kn.  Comparing  rhese  results  with  Table  13  it  is 
obvious  that  the  dynamic  oehavior  of  both  the  linear  and  the 
non-linear  mcdel  are  essentially  the  same. 

Due  to  the  initially  large  plane  deflections  the  speed 
decreases  by  more  than  1  kn  for  all  initially  ordered  speeds 
until  the  ordered  deptn  is  reached.  At  this  time  the  plane 
deflections  get  small  and  the  submarine  accelerates  to  the 
original  speed.  As  predicted,  the  initial  plane  angles  are 
increased  with  decreasing  speed. 

Although  the  reduced  speeds  cause  slightly  increased 
overshoots,  the  designed  controller  is  still  able  to  satisfy 
all  the  design  specifications. 


2.  Turns 


One  cf  the  more  difficult  maneuvers  with  a  submarine 
is  to  change  the  depth  and  the  course  simultaneously.  It  is 
of  interest  to  study  the  behavior  of  the  submarine  in  a  turn 
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first  without  depth  change,  because  the  flow  dynamics  in  a 
turn  change  significantly. 

The  transverse  motion  of  the  submarine  causes  cross  flow 
velocities.  If  the  submarine  had  no  appendices  this 
crossflow  would  produce  no  pitch  or  other  forcing  moment  on 
the  hull.  But  due  to  the  conning  tower  and  other  appendices 
the  flow  over  the  upper  portion  of  the  hull  slows  down.  A 
wake  shed  from  the  conning  tower  also  effects  the  velocity 
distribution  aft. 
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Table  14  -  100  ft  depth  change 
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Figure  109  -  RON  1 
100  FT  DEPTH  CHANGE,  7  KN 
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Figure  110  -  RUN  1 
100  FT  DEPTH  CHANGE,  7  KN 
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Figure  111  -  RON  1 
100  FT  DEPTH  CHANGE,  7  KN 
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Figure  112  -  RUN  2 
100  FT  DEPTH  CHANGE,  9  KN 
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Figure  113  -  RON  2 
100  FT  DEPTH  CHANGE,  9  KN 
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Figure  114  -  RUN  2 
100  FT  DEPTH  CHANGE,  9  KN 
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Figure  115  -  RON  3 
100  FT  DEPTH  CHANGE,  11  KN 
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Figure  116  -  RUN  3 
100  FT  DEPTH  CHANGE,  11  KN 


196 

I 


kiir-H  I,  trif  imoh 


PLANE  DEFLECTIONS  -  10°  PER  INCH 


a 

LT 


TIME  -  80  SEC  PER  INCH 


Figure  117  -  RON  3 
100  FT  DEPTH  CHANGE,  11  KN 
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Figure  118  -  RON  4 
100  FT  DEPTH  CHANGE,  13 
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Figure  119  -  RUN  4 
100  FT  DEPTH  CHANGE,  13  KN 
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Figure  120  -  RON  4 
100  FT  DEPTH  CHANGE,  13  KN 
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Figure  121  -  RUN  5 
100  FT  DEPTH  CHANGE^  15  KN 
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Figure  122  -  RUN  5 
100  FT  DEPTH  CHANGE,  15  KN 
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Figure  123  -  RON  5 
100  FT  DEPTH  CHANGE,  15  KN 
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subaarines  tte  performance  of  the  controller  even  for  hard 
turns  is  very  satisfactory. 

After  both  maneuvers  have  been  studied  independently 
the  combinaticn  of  them  leads  to  Fig.  136-141  and  Table  16. 
It  is  obvious  that  both  parts  of  the  combined  maneuver 
influence  each  other.  The  time  required  to  reach  the  ordered 
depth  increased  with  increasing  rudder  angle.  At  the  same 
time  the  yaw-angle  reached  after  400  s  decreased  by  about 
20°  compared  to  a  turn  without  depth  change.  The  overshoot 
experienced  stayed  still  within  the  design  specifications. 


Table  16  -  Course  and  depth  change  of  100  ft 


Rudder 

angle 

in  ° 

Depth 

reached 

in  sec 

Max 

pitch 

in  ° 

Overshoot 

in  ft 

Max  speed 

deviation 

in  kn 

Yaw 

in  ° 

10 

8S 

1.54 

2.97 

2.08 

70.7 

20 

S4 

2.51 

3.  32 

2.94 

102.8 

30 

98 

3.27 

4.77 

3.9  1 

120.3 

3 •  Cut  cf  trim  conditions 

All  test  runs  performed  indicated  that  the  speed  was 
reduced  due  to  additional  drag  introduced  by  all  maneuvers. 
Frcm  Fig.  1C6  one  must  conclude  that  the  submarine  is  not 
in  trim  during  maneuvers.  This  out  of  trim  condtion  has 
clearly  an  effect  on  the  performance  of  the'  submarine.  But 
as  the  contrcller  is  capable  of  handling  these  deviations  it 
is  not  advisable  to  try  to  trim  the  boat  temporarily  during 
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Figure  124  -  COURSE  CHANGE,  HORIZONTAL  PLANE 
RODDER  ANGLE  =  10°,  SPEED  =  11  KN 
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Figure  125  -  C0DRS2  CHANGE 
rudder  angle  =  IQO,  SPEED  =  11  KN 
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Figure  126  -  COURSE  CHANGE 
RUDDER  ANGLE  =  10°,  SPEED  =  11  KN 
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Figure  127  -  COURSE  CHANGE 
RODDER  ANGLE  =  10°,  SPEED  =  11  KN 
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Figure  128  -  COURSE  CHANGE,  HORIZONTAL  PLANE 
RODDER  ANGLE  =  20°,  SPEED  =  11  KN 
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Figure  129  -  COURSE  CHANGE 
RODDER  ANGLE  =  20^,  SPEED  =  11  KN 
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Figure  130  -  CODRSE  CHANGE 
RUDDER  ANGLE  =  20°,  SPEED  = 
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Figure  131  -  COURSE  CHANGE 
RUDDER  ANGLE  =  20°,  SPEED  =  1 1  KN 
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Figure  132  -  COORSE  CHANGE,  HORIZONTAL  PLANE 
RODDER  ANGLE  =  30°,  SPEED  =  11  KN 
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Figure  133  -  COORSE  CHANGE 
RODDER  ANGLE  =  30°,  SPEED  =  1  1  KN 
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Figure  134  -  COURSE  CHANGE 
RUDDER  ANGLE  =  30®,  SPEED  =  11  KN 
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Figure  135  -  COORSE  CHANGE 
RODDER  ANGLE  =  30®,  SPEED  =  11  KN 
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Figure  136  -  CODRSE  AND  100  FT  DEPTH  CHANGE 
RODDER  ANGLE  =  10®,  SPEED  =  11  KN 
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Figure  137  - 
RODDER 


COORSE  AND  100  FT  DEPTH  CHANGE 
ANGLE  =  10°,  SPEED  =  11  KN 
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Figure  138  - 
RODDER 


COURSE  AND  100  FT  DEPTH  CHANGE 
ANGLE  =  20°,  SPEED  =  11  KN 
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PITCH  -  2“  PER  INCH 


Cl 


TIME  -  80  SEC  PER  INCH 


Figure  139  -  COURSE 
RUDDER  ANGLE  = 


AND  100  FT  DEPTH  CHANGE 
200,  SPEED  =  11  KN 
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DEPTH  -  17  FT  PER  INCH 


na 


TIME  -  80  SEC  PER  INCH 


Figure  140  - 
RUDDER 


COURSE  AND  100  7T  DEPTH  CHANGE 
ANGLE  =  30°,  SPEED  =  11  KN 
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PITCH  -  2®  PER  INCH 


n 

C) 


TIME  -  80  SEC  PER  INCH 


Figure  lU 1  - 
RUDDER 


COURSE  AND  100  FT  DEPTH  CHANGE 
ANGLE  =  300,  SPEED  =  1 1  KN 
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maneuvers . 

To  demonstrate  the  influence  of  the  out  of  trim  conditions 
100  ft  depth  changes  were  simulated.  First,  the  correct 
forward  and  aft  trim  for  the  ordered  speed  of  11  kn  was 
maintained  vihile  the  submarine  was  too  heavy  or  too  light 
overall.  Second,  the  correct  overall  weight  was  used  and 
forward  and  aft  trim  changed  by  the  same  amount,  but 
opposite  in  sign. 


Table  17  -  Incorrect  overall  weight 


Height 

error 

in  lbs 

Depth 

reached 

in  sec 

overshoot 

in  ft 

Steady-state 

error  in 

pitch 

in  ° 

depth 

in  ft 

+  2000 

€4.0 

3.72 

0.035 

0.440 

+  1000 

€5.0 

3.41 

0.018 

0.210 

0 

66.0 

3.  10 

0.000 

0.000 

-1000 

€7.5 

2.76 

-0.018 

-0.315 

-2000 

89.0 

2.45 

-0.037 

-0.516 

Tables  17  and  18  show  the  influence  of  trim  changes.  A 
change  in  the  overall  weight  of  2000  lb  and  a  transfer  of 
800  lb  from  the  aft  to  the  forward  trim  tank  or  vice  versa 
is  required  for  the  adjacent  2  kn  intervals.  The  overshoot 
behaves  as  expected,  i.e.  when  the  boat  is  too  heavy  overall 
or  the  bow  is  too  heavy  the  overshoot  increases.  Vice  versa 
the  overshoot  decreases  when  the  submarine  is  too  light. 
During  all  maneuvers  the  speed  decreases,  which  causes 
according  to  Fig.  108  too  heavy  overall  and  too  heavy 
forward  conditions.  This  explains  the  increased  overshoot 
observed  for  the  depth  changes  simulated  before. 
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Table  18  -  Incorrect  forward  trim 


Forward 

Eepth 

Steady-state  error  in 

trim 

reached 

overshoot 

pitch 

depth 

in  lbs 

iD  sec 

in  it 

in  ® 

in  ft 

+  800 

83.0 

4 . 17 

-0.067 

0.790 

+  400 

85.0 

3.63 

-0.034 

0.380 

0 

86.0 

3.  10 

0.000 

0.000 

-400 

87.0 

2.56 

0.034 

-0.434 

-800 

SO.O 

2.00 

0.067 

-0.840 

A  new  feature  experienced  during  these  trials  is  a 
steady-state  error  in  depth  as  well  as  in  pitch  after  the 
submarine  has  regained  its  original  ordered  speed.  These 
steady-state  errrors  are  small  for  relatively  great  cut  of 
trim  conditions  as  2000  lb  in  overall  weight  and  800  lb 
heavy  forward.  If  the  disturbances  become  greater  it  is 
clearly  necessary  to  prevent  these  steady-state  errors. 
There  are  basically  two  alternatives.  First,  correct  the 
trim.  Second,  increase  the  sensitivity,  i.e.  the  gains  of 
the  compensator  for  small  depth  and  pitch  errors. 
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r  investigation  was  a  multivariable 
with  six  degrees  of  freedom.  Linearizing 
upled  the  vertical  plane,  now  a 
ultivariable  linear  system.  Comparative 
open-loop  linear  and  non-linear  model 
ilization  of  the  linear  model  for  the 
ive  was  to  design  a  depth  controller  with 
of  feedback  states.  Two  approaches  using 
function  matrices  were  investigated. 

ECOUPLING 

wn  that  the  steady-state  decoupling  scheme 
y  applied  to  this  submarine  problem. 

feedback,  steady-state  decoupling  depends 
ntegrators  the  number  of  which  is  easy  to 
lying  the  final  value  theorem.  But  the 
he  compensator  is  a  matter  of  trial  and 
ingle  loop  techniques,  here  the  root-locus 
applied  systematically,  it  still  is  a 
terative  process  which  requires  either 
trials. 


Design  for  stability  and  desired  dominant 
possible,  but  only  simulation  verified  the 
assumption.  As  the  desired  characteristics  fixed  t 
the  compensator  matrix  larger  depth  changes  led  to 
plane  deflections.  Therefore  limiters  were  reguir 
the  planes  within  the  physical  constraints. 


roots  was 
dominance 
he  gain  of 
excessive 
ed  tc  keep 
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The  ccntrcller,  designed  on  the  basis  of  the  linear 
model  dynamics  for  9  kn,  showed  the  following  performance, 
when  implemented  in  the  non-linear  model: 

a. )  The  design  specifications  for  a  speed  range  from 

7  -  15  kn  were  satisfied. 

b. )  Out  cf  trim  conditions  led  to  small  steady-state 

errors . 

c. )  When  operated  in  depth  keeping  mode  in  a  turn  the 

ccntrcller  performed  satisfactorily. 

The  compensated  system  is  rather  insensitive  to  speed. 
Thus  all  problems  related  to  gain  switching,  like  chattering 
and  discontinuities  in  command  plane  angles,  are  avoided. 
This  is  especially  important  because  the  speed  changes 
significantly  during  maneuvers. 

The  design  was  only  concerned  with  the  derivation  cf  the 
transfer  functions  for  the  cascaded  compensator.  When 
implemented  in  hardware  it  has  the  following  desirable 
features: 

✓ 

*  Minimal  instrumentation 

As  no  rate  information  is  required  only  simple  sensors  for 
depth  and  pitch  angle  are  necessary,  i.e.  no  inertial 
guidance  system. 

*  Low  cost,  weight,  and  size 

The  simplicity  of  the  compensator  transfer  functions  makes 
them  easily  realizable  in  physical  hardware  at  low 
manufacturing  cost.  Weight  and  size  requirements  are  very 
small,  another  important  factor  especially  for  small  coastal 
submarines. 

*  Reliability 

The  automatic  controller  can  be  realized  with  a  set  of 
physical  components  with  a  well  known  high  reliability.  High 
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ccnponent  reliability  and  a  small  number  of  components  will 
generally  result  in  a  high  system  reliability. 


E.  TOTAL  CECODPLING 


As  for  the  steady-state  decoupling  scheme  the  tctal 
decoupling  is  based  on  the  linear  model.  The  determination 
of  the  compensator  matrix  is  an  easy  mathematical  process, 
once  the  desired  response  dynamics  have  been  defined  in  form 
of  an  overall  transfer  function  matrix.  At  the  first  glance 
this  method  seems  to  be  very  appealing  but  it  has  some 
limitations : 
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with  speed,  the  compensator  theoretically  covers  only  a 
specific  Sfeed.  Test  runs  showed  that  a  compensator 
calculated  fcr  one  speed  was  capable  of  covering  a  very 
small  speed  interval.  In  order  to  cover  a  large  speed  range 
as  required  in  our  problem  the  compensator  has  to  be  very 
complex  because  the  gain  and  the  root-locations  have  to  be 


adjusted  in  very  small  step  sizes.  The  influence 
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this  investigation  it  is  c 
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cascade  compensation,  using  the  steady-state 
method,  is  possible  and  practical  for  automatic 
depth  control  cf  small  submarines. 


decoupling 
pitch  and 


229 


'M'  . 


i 


APPENDIX  A 


A.  AXIAL 

m(u  - 


STANDARD  EQUATIONS  OF  MOTION 


FORCE 


+  %  I 


^  fx  .  '  U 

L  u 


^  X  .  ■  U  +  X  ^  wq  ’ 


] 


*  i‘’  [".u 


• 


V*  +  X  *  w^l 

VV  ww  J 


^  I  •  6,  '+  X  -  4.'  +  X 


4r6r  ''r  ^  ''  4s68  ®8  4b4b  *b 


Mn'  -  1)V 


'I"’  -  ')  V 


E\V.  sin  9 
1 
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.A 


•» 


fi 


B.  LATEBAL  FOECE 
m(v  -  wp  +  ur)  -  "2  ^ 

♦f‘‘ [%,  ■‘>'>*’fplpl  'pIpI] 

[v  .  ■  V  *  Y  ^  •  wp  t  Y,|,|.  V|  |(v-  t  w‘)i  I  |r|  ] 

[y^  •urtY|^|,^  'ulrUrtY^  'up] 

*i‘‘  rn'  '  <'>'  -  '>“■■ 

[y  .  '  u’  *  Y  ^  '  uv  t  Y  ^1^1  '  v|(v‘  * 

'  ‘r 

'In'  -1)5, 

+  |/'  Y  '  (n'  -  1)  uv 

■^f^^^v|v|n'  ’  (n'  -  1)  v|(v^  +  1 

+  ■?  Y  '  wv# 

C  v/V 

+  *2  (^y)  (-W')  Sin.cot  f  Multiplied  by 

^  •  u 

U 

I  for  large  angle 

+  EWi  sin(;)  cos  0  attack  near  -90 
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NCBHAL  FCBCE 


m(w  -  uq  +  vp)  =  7  Z  .  '  q 


l‘*  1 

•  *  +  Z  ' 

rp 

rp#] 

NoC«  1 

0  .  1 

fz  . 

•  w  +  Z  '  vr#  +  Z  • 

vp  +  az  ' 

L  w 

vr 

vp 

•  vp 

|z* 

[s 

1,U. 

*  u  1  ql  6 s  + 

^wui?r' 

^  1"' 

*  (n*  -  1)  uq 

i‘‘\ 

*  +  Z  *  uw 

w 

^  ^  wjwl 

’  wl(v*  +  W^)‘ 

*  ■i‘‘  I'"''’"*  ”•■’'‘1  *  ^v' 

] 

*i‘‘  w„.  'In'-  l)uw  +  Z^|„|„, 


V  Sin  cot 


+  SW.  cos  9  cos  ^ 


H  Multiplied  by 

u 

U 

for  l.^rge  angles  of 
attack  near  -90^ 

r^ote  1 

when  not  multiplied  by 
u  add  to  Z  * 

D 


232 


\ 


D 


ROLLING  i!CMENT 


.  (I^  -  ly)  <ir 


0  ,6 
7  ^ 

[k  .  ' 
L  p 

'  p  +  K  ' 

^  qr 

'  pIpI] 

f‘‘ 

up  +  K  ^ 

•  ur  +  K  .  •  V  +  K  ' 

V  wp 

wpj 

[•<.• 

u'  +  K  ' 

V 

“^■•■^vlvl  "''ll''*  +wM^|] 

|x“ 

K 

vw 

•  vw 

I'’ 

u^K  . 

6r 

’  ^ 

• 

+  Bs  g  flin  C09  0 
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E.  PITCHING  MOMENT 


Note  1 

V  ^  [‘‘’‘  4  '  ^  rr  '  ,p  '  *-P  +  AM  '  rp#] 

*  [^  q  '  “'I  +  ^  |q|6,  '  «lql4s  +  M  1^1^  •  |fv'  +  w*^^lq] 

+  -|z*  [m  ^  '  w  +  M  '  vr#  +  M  ’vp#] 

+  l"*-'^qa'  '  <"'  - 

+  1^=*  [m  ,  '  u^+  M^  '  UW+  M  •  w|(v>+w»)^|] 

^1"'  [^6s  '^■^'^‘ob  '^] 

+  2^"*  wn*  ' 

+  f  w|wla'  '  <"'  -  1)  w[(v2  +  ^2)^  I 


"  I'"  '("’* 


+  Bz  ^  sin  0 


-  EW  ^  x^.  cos  0  cos  ^ 


?  Multiply  by 
large  angles  o 
near  -90'^ 


Note  I 

when  not  multiplied 
by  ^  add  to  M  ' 

U  ^P 
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F.  YAWING  MCMZNT 


.  t.t.  -I  IP,  [«f  >■] 


2  y 


l‘‘  1 

'  ur  +  N 

1*'  i 

N 

P 

•  up  +  N  , 

-f'' 

N  , 
rn* 

'  {o’  -  1 

0  a  3 

2 

'  u*  +  N  , 

D  a3 

2  ^ 

u*  N 

6r  ’‘r 

u*  N 

6rn'  ’ 

N 

vn 

,  '  (n*  - 

Nvlvl  '  ^  0 


+  0/^  N  ,  ,  ,  •  (r^'  -  1)  vl(v'  +  wM^ 

2  vlvj  o' 


+  ^  N  '  wv/f 

^2  wv 


+  L'>'''-  X.-  cos  0  sin  '^) 
»  ti 


f  Multiply  by  -  for  large 
angles  of  attack  near  -90* 
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G.  AOXILIARY  EQUATIONS 


({>  5  p  t  sin0 

0  =  (q  -  ^Jco  5  9  5  in  (})  )  /  CO  5  (J) 

^  =  {  r  +  0  sin  4>)  /  cos  0  CO  3  4) 

Xo  =  ucoaOcosv^  +  v  (sin^)  sin9  cos^  -  cos^  sin  ij)) 
+  w  (sincf)  sin^  +  coscj)  sin©  cos 


Yo  =  u  cos  9  sin  ip  +  v  (cos({)cos;{»  +  sin<})sin0  sin  ) 
•f  v/(cos(J)sin0  sinijj  -  sin  ^  co  s  ip  ) 

Zo  =  -a  sin9  +  vcos9  sincj)  +  W’cos0cos4) 

U  =  (u’  + 


(f'x)p 


f.'  u‘| 

[  ^i'  +  ig'  n*  +  ag*  n‘^  J 

[v  +  b3'n-+b3-  n'^] 

f/U'l 

[  c;  +  c,'n'  +  C3-n-'] 

[  d,'  +  d,-n-  +  d3-n-^] 

when  <  n* 

when  kg  <  n*  <  kj 


when  kg  <  n*  <  kg 


when  n'  <  kg 


•  t  , 

V  I  t 

^2i  ^3 

•  *  i 

Cli  C2,  C3 

I  I  • 

di,  do,  dj 


Sets  of  non-d l'’.er>s lor.a I  cocrflclents  'iscd  In  the  pro¬ 
pulsion  oq^iatlon  above,  "'h^'  set  vhlch  will  be  In  etTect 
at  any  tl.T.c  d.irir.,’  a  slnilatci  r^ancuver  vilL  depend  on  the 
value  of  n  and  tl*.e  nur/ners 
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APPENDIX  B 


NOMENCLATURE 


SYMBOL 


DEFINITION 


,  A  dot  over  any  symbol  signi¬ 

fies  differentiation  with  res¬ 
pect  to  time. 

B  Buoyancy  force  which  Is  posi¬ 

tive  upwa r ds  . 


n 


Si 

u 


U 


V 


V 


Mass  of  the  submarine  Including 
the  water  In  the  free  f3oorl<ng 
spaces • 

Overall  length  of  the  submarine 

Linear  velocity  of  origin  of 
body  axes  relative  to  an  earth- 
fixed  axis  system. 

Component  of  U  along  the  body 
x-axls . 

Component  of  U  along  the  body 

y-a  x Is . 

« 

Component  of  U  along  the  body 

z- ax 1 s . 
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Comsiand  speed:  A  steady  v.ilue 

of  u  for  a  given  propeller  rpm 
whenc^»^  and  control  surface 
angles  are  zero*  Sign  changes 
with  propeller  reversal. 

Longitudinal  axis  of  the  body 
fixed  coordinate  axis  systeta. 


Transverse  axis  of  the  body 
fixed  coordinate  axis  systeo. 


s 


X 


o 


Vertical  axis  of  the  body  fixed 
coordinate  axis  systeci. 

Distance  along  the  axis  of  an 

earth-fixed  axis  aysten. 


Distance  along  the  y^  axis  of  an 
earth-fixed  axis  systca. 

Distance  along  the  axis  of  an 

earth-fixed  axis  system. 


Component  of  angular  velocity 
about  the  body  fixed  x-axls. 


Component  of  angular  velocity 
about  the  body  fixed  y-axis. 


Component  of  angular  velocity 
about  the  body  fixed  z-axls. 


z 


E 


The  z ^coordinate  of  the  center 
of  buoyance  (CB)  of  the  subma¬ 
rine. 
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1 


o< 

An  Ic  of  attack. 

An$:le  of  drift. 

Sb/  Dfa 

Deflection  of  bowplane  (  or 

a  a  1  Ip  lane  ) 

Sr  ;  Dr- 

Deflection  of  rudder. 

Ss  j  Ds 

Deflection  of  aternplanc. 

f 

O 

0 

$ 

</> 

C 

The  ratio  u^/u. 

PI t ch  angle . 

Yaw  angle. 

Rol 1  angle. 

Maas  density  of  sea  water. 

«i 

Weight  of  water  blown  from  a 

particular  ballast  tank  Idcnt- 
Iflcd  by  the  integer  assigned 

to  the  Index  1. 

6J 

Angular  velocity. 

t 

Tine  . 

^tl 

Location  along  the  body  x-axis 

of  the  center  of  nass  of  thP  1 

ballast  tank  when  this  tank  is 

filled  with  sea  water. 
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■f 

I'  f 


,  4  '' 


i 


t 


I 


X 


Propulsion  force  (see  auxi¬ 
liary  equations  and  relation¬ 
ships)  . 

Mooicnt  of  Inertia  of  a  sub¬ 
marine  about  the  x-axla. 


I 


7 


Moaent  of  Inertia  of  a  sub¬ 
marine  about  the  y-axls. 


I 


z 


Moment  of  inertia  of  a  sub¬ 
marine  about  the  r-axis. 


s'*  S|p|'>Sr' 

S'  S*  s'.  Sp'.  K/. 

^vjvj'.  <vw'>  i^Sr' 


Non- d ine ns i on j 1  constants  each 
of  which  is  assigned  to  a  parti¬ 
cular  force  tern  In  the  equation 
of  motion  about  the  body  x-axis. 


M. 

q 


H  •  , 

r  r 


AM 


rp 


M. 

I 


q  I  43  * 


Ml  I  ',H  ',M  ’.M**, 

|w|q  V  vr*  VO*  nn*  ** 


vr  vp  qn 


M  *  ,  M  I  1  ,  Mt  I  •  ,  M  ’  ,  H  ’  ,  M  '  , 
V  w  w  w  • ww  *  vv  63 


M..I. M,  .* 


Non-dimensional  constants  each 
of  which  Is  assigned  to  a  parti¬ 
cular  force  tern  In  the  equation 
of  motion  about  the  body  y-axls. 
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%<,’•  "p’-  "lr|4r’*  ^|v|r’ 


N  «  vj  •  W  *  S’  *  S  *  S  * 
p  •  •  ^wp  '  ^rn'  •  ’  V 


V  V 


’•  ^Sr’*  ^rn*'*  •'''vn*’*  ^|v!n-’' 


Son-d Incns tona 1  constants  each 
of  which  Is  assigned  to  a  parti¬ 
cular  force  tern  In  the  equation 
of  notion  about  the  body  r-axls. 


5  • 

vv 


X  * 

<1<1 


X  ’ 
r  r 


rp 


X.  ’ 
u 


X  • 

V  r 


X  * 
wq 


Y*Y*Y*Y  * 

*uu  •  ^vv  •  ^ww  •  ^6  r6  r  ’  5  s6  a  * 
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APPENDIX  D 


DSL-PROGRAMS  AND  BLOCK  DIAGRAMS 


A.  COMPENSAIED  LINEAR  MODEL 
(STEADY-STATE  CECOUPXING) 


B.  COMPENSATED  LINEAR  MODEL 
(TOTAL  DECOUPLING) 


C.  COMPENSATED  NON-LINEAR  MODEL 
(STEAEY-ST ATE  DECOUPLING) 
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